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EQUILIBRATING EFFECT OF MAXWELL-TYPE BOUNDARY 
CONDITION IN HIGHLY RAREFIED GAS 


HUNG-WEN KUO 


Abstract. We study the equilibrating effects of the boundary and inter- 
molecular collision in the kinetic theory for rarefied gases. We consider the 
Maxwell-type boundary condition, which has weaker equilibrating effect than 
the commonly studied diffuse reflection boundary condition. The gas region 
is the spherical domain in d = 1,2. First, without the equilibrating ef¬ 
fect of the collision, we obtain the algebraic convergence rates to the steady 
state of free molecular flow with variable boundary temperature. The conver¬ 
gence behavior has intricate dependence on the accommodation coefficient of 
the Maxwell-type boundary condition. Then we couple the boundary effect 
with the intermolecular collision and study their interaction. We are able to 
construct the steady state solutions of the full Boltzmann equation for large 
Knudsen numbers and small boundary temperature variation. We also estab¬ 
lish the nonlinear stability with exponential rate of the stationary Boltzmann 
solutions. Our analysis is based on the explicit formulations of the boundary 
condition for symmetric domains. 


1. Introduction 


In kinetic theory, a fundamental and central issue is the equilibrating effects of 
the boundary and intermolecular collision. In the present study we consider the 
Maxwell-type boundary condition: 


( 1 ) 


/ Ott \ 2 

Fiy,C,t) = «(y) jF{y,t)MT(y)(C) 

+(1 - a(y))^(y, C - 2(C • n)n, t), y S C • n > 0, 


jF{y,t) 


= j -C:,-nF{y,C^,t)dC^ 


boundary flux of E, 


C.-n<0 


where F is the velocity distribution function of the gas particles, C is the microscopic 
velocity, T(y) is the boundary temperature at the boundary point y, n is the unit 
normal vector at the boundary, pointing to the gas region U, a(y) (0 < Q;(y) < 1) 
is the accommodation coefficient, and Mt is the Maxwell distribution: 

g 2RT 

FIt{C) = -Ti F : Boltzmann constant. 

{2Tr RTy 

The case a = 0 is called the specular reflection boundary condition, which has 
no equilibrating effect. The case a = I is called the diffuse refection boundary 
condition, which has strong, direct equilibrating effect of the boundary thermal 
information on the gas flows. In this paper we assume that the accommodation 
coefficient a, 0 < a < I, is constant. Our purpose is to study the equilibrating 
effect of the Maxwell-type boundary condition, the dependence of the process of 
convergence to steady states on the accommodation coefficient a, 0 < a < 1. The 
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equilibrating process also depends on the geometry of the boundary. Our analysis 
demands the quantitative method for the study of particle propagation. For this, 
we will focus on spherical symmetric domains: 

Z? = {x e : |x| < 1} , 

in space dimension d = 1,2. This allows us to use the stochastic formulation of our 
previous works |10j and m , which provides an explicit description of the evolution 
of the free molecular flow. 

We decompose the microscopic velocity C = (Cij C 2 , Cs) £ into 

(2) 4 = (Ci,...,Cd)eR", »7 = (Cd+i,..-,C3)GR'-", 

and rewrite the Maxwell-type boundary condition © as 


(3) 


/ 27r \ 2 

Fiy,C,t) jFiy,t)MT(y){C) 

+(1 - a)F{y, I - 2(1 • n)n, -q, t), y G dD,^- n> 0, 
jpiy, t) = J -I, • nJ^(y, C, t)dC* : boundary flux of F, 


€.n<0 

To focus on the equilibrating effect of boundary, we first consider the free molec¬ 
ular flow: 

' X! ^ X G C G t > 0, 


(4) 


dt 


5(x, C, 0) = C), X e D, C G R^, 

/ 27r \ ^ 

g{y,C,t) = a jgiy,t)MT{y)iC) 

+(1 - a)g{y, | - 2(| • n)n, 77 , t), y G 5£), | • n > 0, 

Js(y> t) = / -I* ■ ng{y, C*, t)dC* : boundary flux of g. 

«.n<0 

The equation for the steady state of the free molecular flow is: 
dS 

5 = 5(x,C), XG A CGR^ 


(5) 


/ S{x, (^)dxd(^ = 1 : unit density, 
\F\ JdxR3 

/ 27r \ ^ 

5(y,C) = a:(^-^^^j js{y,t)MT(y){C) 


+ (1 - Q!)S'(y, I - 2(1 • n)n, q), y G dD, ^ ■ n> 0, 
js{y) = / -C ■ nS'(y, C*)c?C, : boundary flux of S. 
Here we take general initial data gin with finite weighted norm: 


( 6 ) 


5i„(x,C) G L'^f, M > 4, 


|| 5 „||^oo^ = ||g„||^ = esssup (1 +|C|)''|ffm(x,C)|, 


where the choice of /r > 4 implies that 


f 

-I • n 

/^■n<0 

(1 + ICI)^ 


dC < 00 . 
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The boundary temperature variation is assumed to be bounded: 

0 < T, = inf T(y) < T* = supr(y) < oo. 

dD gjj 

The diffuse reflection boundary condition has strong and direct equilibrating ef¬ 
fect, and as a consequence, the convergence to steady state is of the rate of d 
the space dimension, nniiii]- The Maxwell-types boundary condition yields even¬ 
tually the same rate, with intricate dependence on the accommodation coefficient 
a. 

The following Theorem shows the convergence to the steady state of free molec¬ 
ular flow S. 


Theorem 1.1 (Main Theorem for Free Molecular Flow). For G and 

/i > 4, the solution of o satisfies 


g{x, C, t) - p^S{x, C) = 0(1) hinW 


M(C) ( 1 -a)^' 


°°’>^\y{l + aty ^ (1-blCI)'^ j 
1 


-f 


(1_P|^|)m^{I€I<P^}/’ 


P*^ — 


\D\ 


C) dxd^ , 


’Dxl 


for any small e, 0 < e < 

Theorem ll.il immediately implies the following convergence of g: 

Corollary 1.2. For any small e, 0 < e < g converges to pt,S in for 
1 < P < oo: 


\\g{x,C,t) - P*S{x,C)\\lI^ = 0(1) ( 

Consequently, there exists Ca,e > 0 such that 


1 


{at + l)'^ 


-f(l-a)*02+ 


1 


{t + 1) 


(1-0^)' 


\\g{x,<!^,t) P:i,S{x, <^)\\ Ca,e \\gin\\aci,ij, ( d 


In particular, the coefficient Ca,e = 0(1) when a = 1. Hence, we may let e ^ 0 to 
obtain the optimal rate for diffuse reflection boundary condition: 


\\g{x,C,t) - p,S'(a;,C)||LP = 0(1) ^ • 

V(^ + 1)V 

After studying the boundary effect of Maxwell-type condition for free molecular 
flow, we continue study the additional equilibrating effect of the collision in rarefied 
gas flow. We use the Boltzmann equation to model gas with intermolecular collision. 
Consider the initial-boundary value problem of the Boltzmann equation: 


( 7 ) 


" 3F ^ BF 1 

^ + xeiicR^ CeK3, t >0, 

7^(x,C,0) = Fi„(x,C), X e D C M"*, c e 
Maxwell-type boundary condition ©■ 
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where k is the Knudsen number which measures how rarefied the gas is, and Q{-, •) 
is the collision operator, a symmetric bilinear operator. 


Qig, hKC) = \ j {giC'MC) + HOgiO - giCMC) - HOgiO) 

S2xR3 

xB{9,\c-C\)dndC, 


where 


fc' = c-( 

(C-C*)-f^)^^ 

)«;=<.+ 

{ic-a-^y 


cosd = 


c-c 


K-C*l 


and B is the collision kernel which is determined by the interaction potential be¬ 
tween two colliding particles. Throughout this paper we assume an inverse power 
hard potential with Grad’s angular cut-off or hard sphere. Under this model. 
Bid, 1C* - Cl) ^ 1C - C*l^l cos6>|, for some m > 4. 


Maxwell molecule 

hard potential 

hard sphere 

u = 4 

4 < u < 00 

u = 00 


By nondimensionalization, m, we may assume, without loss of generality, that 
0 < T* < T* = 1 and the total density unity: 

(8) / F,„(x,C)dx(iC = 1- 


For convenience, denote the Maxwellian Mt*(C) = Tfi(C) = (tt) 2 exp(—|CP) 
simply by M(C). 

Conventionally, to linearize the Boltzmann equation, we expand F around M, 
F = M + '/Mf. The resulting equation for the perturbation / is 


dt 




dxi 


K kVM 


Q 



where the linearized collision operator L is defined as 


Lf 




Q VM, 


,m) 


and the linearized Boltzmann equation is 


dt 




dxi 


1 


-L/ = 0. 


For the intermolecular force model we consider, inverse power hard potential with 
Grad’s angular cut-off or hard spheres, L can be decomposed as the difference of 
an integral operator K and a multiplicative operator v. 

L = K-v, [Kf){C) = J^^KiC,C)fiC)dC., (^^/)(C) = ^^(C)/(C). 

u(C) is a positive function with non-zero infimum: 


uq = inf u(C) > 0. 

CGR3 

However, as M does not satisfy the boundary condition Q, this linearization is not 
natural for the present situation, where the boundary effect is significant. Instead, 
we expand around the stationary free molecular flow S under the Maxwell-type 
boundary condition and write F = S + VMf- Because S does not satisfy the 
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Boltzmann equation some more source terms are introduced to the equation 
for / : 


(9) 


dt 





+ 


dxi 




M + y/Mf, S-M + 



On the other hand, we have the important property that S satisfies the Maxwell- 
type boundary condition. As the boundary condition is linear and homogeneous, 
the boundary condition for y/Mf remains the same: 


( 10 ) f{y,C,t)VWO 

MT(y)iC) 

+ (1 - a)/(y, I - 2(1 • n)n, rj, t)\jM[C) 

y e dD, I • n > 0. 

Therefore, we will consider the initial-boundary value problem of equation ([9]) with 
boundary condition (US and initial data: 


/zn(x,C) 


f^n(x,C) - 5'(X,C) 


There is a trade-off between more complicated boundary condition and extract 
interior source terms. We choose the latter since the effect of Maxwell-type bound¬ 
ary condition has been well analyzed for the free molecular flow. Moreover, this 
linearization is physically natural for large Knudsen number. The extra interior 
source terms can be handled by the standard iteration scheme. Therefore, we will 
first consider m as our linearized equation: 


( 11 ) 


dt 




Maxwell-type boundary condition (uni). 


/DxR3 


VA^/„(x,C)dxdC = 0. 


From F = S + \fMf, the zero total mass condition J y/MfindxdC = 0 is a con¬ 
sequence of nondimensionalization, (jS]), rather than an additional constraint. To 
study the problem (ED, we make two reductions: first we reduce (9t-l-X) Cidxi 
to {dt -f X) Cidxi + ^), and then reduce {dt + CAi + to {dt + Odxi)- 

For the linear problem ED, we prove an exponential decay to zero of /, Theorem 
o We then use this result to: (i) construct the steady state solution of the 
Boltzmann equation Q , Theorem 11.51 (ii) obtain an exponential convergence to 
the steady state, Theorem ll.61 

Theorem 1.3 (Stability for Linear Boltzmann Equation). Suppose that fin G 
for some constant 7, 0 < 7 < 1, / fm'/MdxdC = 0, and that, for each 
0 < < vq, there exists a positive constant C\ such that for all Knudsen number 
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K with 


( 12 ) 


K> Cl 




logK 


>Ci 




(2|logKuo-u()]|)^°° \ 
a{vQ - u()|log(l - 

_ (2|log[a(uo-u()]|)4°° \ 

u()3/2 a(uo - u()| log(l - a)|4ooy 


ford= 1, 
for d = 2, 


the solution f of (HH) decays to zero exponentially in time: 


(13) 


||/(•,•,^)IL,_^<C||/,„||^_,e- 


for some constant C independent of k, a, T* and . 

Remark 1.4. Our main interest is in the highly rarefied gas, i.e. the case of large 
Knudsen number k. Thus we first study the free molecular flow for collisionless 
gas. Then we consider a perturbation around the steady solution of free molecular 
flow for the Boltzmann equation. To obtain the exponential stability for linearized 
Boltzmann equation, we start with the estimate of free molecular flow which is 
a limiting case of k = oo. Recall that the pointwise estimates of free molecular 
flow depend on the accommodation coefficient a. Theorem 11.11 Consequently the 
magnitude of k is related to that of a, m- 

In this paper, we consider the situation of variable boundary temperature. It 
is a highly non-trivial problem to study the existence of steady solution for the 
Boltzmann equation when the boundary temperature varies. In the case of diffuse 
reflection with variable temperature, the existence of the steady solution in a convex 
domain with dimension less or equal than three was proved by Guiraud [SJ [5] for 
arbitrary fixed Knudsen number. The same result was proved for arbitrary domain 
together with the exponential stability [7]. Moreover, a large data existence result 
was proved by Arkeryd and Nouri for prescribed initial data [3]. In our work, 
with the exponential convergence of linearized Boltzmann equation, Theorem 11.31 
we are able to construct the steady solution as a consequence of time asymptotic 
analysis. To handle the nonlinear term, we require the small variation of boundary 
temperature 1 — T, ^ 1. Then we prove the existence of the steady solution for 
the Boltzmann equation ©■ However, the aim of this paper is to establish the 
equilibrating effect of boundary and collision. The method we used in this paper 
can be seen as an alternative approach to the existence problem of steady solution. 
The quantitative structure of the steady solution on the temperature variation is 
remained for our future research work. 


Theorem 1.5 (Existence of Steady Solution for Full Boltzmann Equation). As¬ 
sume that 1 — T, <C 1 and k ^ 1 satisfies CH). Then the steady state solution $ 
of © exists and satisfies 


(14) 

(15) 


m 


= o(i-r,), 


/UxR3 


^{x, C) \/M{C)dxdC = 0. 


We have already obtained the steady state solution Fao = S' + for full 

Boltzmann equation ©. Moreover, from (ITK1) . f FaodxdC = 1. For general initial 
boundary value problem, we expand F around Fac,: F = Foo + '/Mif. The equation 
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for Ip is 


(16) 


0^ ^ dP; 1 

1 — 1 

V'(x,C,0) = i/’m(x,C) = e 

Maxwell-type boundary condition m 


/DxR3 


■*/'m(x,C)\/M(C)dxdC = 0. 


To reiterate, zero initial total molecular number J ipin'/MdxdC = 0 is a conse¬ 
quence of nondimensionalization (|S]), rather than an additional constraint. 

To show the stability of steady solution Fao , it is sufficient to show that tp decays 
to zero. We establish the exponential decay rate by using Picard iteration with 
the estimate (0 of linearized Boltzmann equation. The following theorem shows 
that the steady solution of the Boltzmann equation is exponential stable when the 
Knudsen number k is sufficiently large, m- 


Theorem 1.6 (Stability of Steady Solution for Full Boltzmann Equation). Suppose 
that (fl^ and 1- Then, for any fixed 0 < ui < uq, the solution ip 

of (fTC)) exists and satisfies 

for a positive constant independent C of a, k, T* and ui. 


The present study focuses first on the equilibrating effect of the boundary condi¬ 
tion by considering the free molecular flows. This approach has been considered for 
the diffuse reflection boundary condition, initiated by m for one space dimension 
and then generalized to higher space dimensions by m for constant boundary tem¬ 
perature, and by m for variable boundary temperature. It also has been studied 
for half space under gravitational force [12]. This approach makes it possible to 
study of the full Boltzmann equation when the Knudsen number is large. Other 
studies consider the case when the Knudsen number is of order one. In other words, 
they consider the case when the collision plays a role at least as important as the 
boundary condition in the equilibration of the gases. Also, other studies consider 
the diffuse reflection boundary condition. For this, see laiiiiii] when the bound¬ 
ary temperature is constant, and miiisiiaiT] for variable boundary temperature. 
For the specular reflection condition, the equilibrating effect has to come from the 
collision. There have been substantial progresses in this regard on the level of the 
Boltzmann equation, see HEKIS], and references therein. 


2. Free Molecular Flow 


2.1. Preliminaries and main results. The steady state solution S of free molec¬ 
ular flow under the Maxwell- type boundary condition m has been constructed 
explicitly, m- 


( 17 ) 


1 °° . / 

^(x,^ ( 

i—1 


27r 


Vi?r(x(i)) 


^T(x(,))(C) 




( jir(x,.|.) j W(x,„.>K.)fcdC 
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where the boundary point obtained by tracing back from the given interior point x 
along the direction — ily: 

^x, ill ^ = X — I X sup {s > 0 : X — G I?, for all s' G (0, s)} , 

X(i) = Yb (x, ify) , 

=|-2(|-n(x(i)))n(x(i)), 

^(fe-l-l) Vb (^{k)-i I 

^fe+i^^fe_2(^fc.n(x(,+i)))n(x(fc+i)). 

Since the domain D is symmetric, we have the following lemma: 

Lemma 2.1. For {x,^) G H x and each k>l, 

|®(/c+l) ^(/c) I “ |^(/c+2) ®(fc+l) I ^ I® 

From the explicit expression (II3, 

(18) ^(x, C) - M(C) = 0(1 - T*)M(C). 

We note that S has constant boundary flux l/Cs- 

[ -^■nS{yX)dC 

J |-n<0 

1 /* ^ f 27r \ ^ 

=—a^(l - (47r)W ■ nM{C)dC = —. 

Os J^iKO Os 

Note that both the evolutionary equation (j?]) and the boundary condition 
conserve molecular number, therefore the total molecular number f ( 7 (x, t)dx(i^ 
is a constant of time. We define the average total density as: 

P* = T^I 9in{y^X)dxdC =-r^ I g(x,C,t)dxdC, 

l-^l JdxR^ \^\ JdxR3 

a constant associated with gin- Due to the equilibrating effect of the Maxwell-type 
boundary condition, one can expect the solution g to approach the steady state 
p*S'. Namely, we expect the function g — p^S{x.X) to decay to zero. Moreover, 
g — p^S satisfies the same evolutionary equation (U) and the boundary condition 
@. Since the space dimension is d and d < 3, it is natural to integrate out the 
extra microscopic velocity degrees of freedom: 


(20) 

g{xX,t) = 

/ (5(x,C,t)-/3*-S'(x,C))di7, 



jR3-d 

(21) 

ff*n(x,|) = 

/ {9^ni^,C) - P*S{K,C))dr], 



yRS-'i 

(22) 

s(x,|) = 

[ S'(x,C)dJ7- 

jR^-d 
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Recall that rj is the last 3 — d components of C, ©■ Since S has constant boundary 
flux, the corresponding boundary flux becomes jg — pt^jCs'- 


j -I ■ ng(x,|,t)d| 


^•n<0 


(23) 


I -l-nf J (5(y,C,i) - P*<5'(y,C))di7 I 


^•n<0 


J -I ■ rig{y,C,t)dC- p*/Cs = jgiy,t) - ^ =jg{y,t)-js, 


in<0 

and the total molecular number becomes zero: 


(24) 


J - p*S{x,C))dxdC = J g(x,|,t)dxd| = 0. 


Dxl 




Moreover, the new functions g(x,t),j(y,t) satisfy equations similar to that for 
the original functions: 


(25a) 


(25b) 


If+ = g = ff(x,et), xGi4cR^ I gM^ t >0, 

XXX,o) = g„(x,|), 

/ 27r \ ^ 

giyX.t) = y jiT{y) ) diy.t)MT(y)X) 

+(1 - a)g{y, ^ - 2(^ • n)n, t), y G dD, | • n > 0, 

/ fi 2RT 

MtX) = / MT{C)dii = - 

J ^ (2ttRTY 

*3-<i ' > 


but with the additional zero total molecular number condition: 

(26) J g{x,^,t)dxd^ = 0, t>0. 

DxR"^ 

Note that Mt{C), the Maxwellian, and MtX), the reduced Maxwellian, are gen¬ 
erally different as functions. To avoid confusion, we always refer to M as the 
abbreviation of MX), not MX). 

For X G D and ^ G R'^, we define Tf, = rb(x, ^) the backward exit time: 

(27) Tt,(x, ^) = sup {s > 0 : X — G D, for all s' G (0, s)} , 
and 


ti = Tt,(x,|) = 


X — X 


( 1)1 


III ’ 


^ __A- ^fe^ _ l^(fc) -^(fe+l)l 

4+1 -r6(x(fe),C )- --. 


From Lemma \TT[ we have tk = 4 for all fc > 2. 
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Suppose that the boundary flux j is given. Then the solutions of the transport 
equation (12511 has explicit form by the characteristic method: 


(28) g{xX,t)- 


(29) 

5(x,|,t) = < 


p,S'(x,C) 

' m-l 

a ^(1 - - kt2) - 

k^O 

_ ^ +((1 - a)™g^n(x(„) -|'”(t-ti - (to - l)t2),^™, J?) 

OO 

-a ^ (1 - af for n < t, 

k—m 

gin{x- ^t,C) - P*S{x,C) iort<n, 

m— 1 

a ^ (1 - afj (x(fc+i) kt 2 ) (l'") 

+(1 - a)™g,„(x(„) - |™(t - ti - (to - 1 )^ 2 ), o for n < t, 

Qin (x-|t,|) fort<rb, 


where 

(30) 


m = 


llitdiiwdi + i 

|X(1) - X(2)| 


and for simplicity of notation we set 

f 2 tt 

The following are our main theorems for free molecular flow, which will be proven 
in the following four sections. 


Theorem 2.2 (Global Existence for Boundary Flux). The solution of (HJ and 
(HSl), with initial data gtn € , exists and is unique for g, > A. Moreover, there 

exists C > 0 such that 

(31) j(y, t) = 0(1) WginWao,), e^^^\ 

where C depends on T* and T,. 


Theorem 2.3 (Decay Rate for Boundary Flux). Suppose that g^ € for some 
constant g> A. Then the boundary flux j{y,t), (l^lTll . satisfies 

+ (1 - «)■*) 

for some constant C depending only on g, T, and T*. 

From (l30l) . for 0 < e < 1 we have m>t^ for |^| > Therefore, Theorem 12.31 
together with (l28l) yield immediately the pointwise convergence of the free molecular 
flow g, Theorem ll.il 


Remark 2.4. When 0 = 1, the case of diffuse reflection boundary condition, 
we have the convergence rate (1 + t)“‘^ of the boundary flux, (TUI |TT]. Roughly 
speaking, the equilibrating effect is mainly from sufficiently many collisions with 
the boundary of diffuse reflection condition when t is large. For Maxwell-type 
boundary condition, we have two possibilities after each collision with the boundary: 
one is diffuse reflection and another is specular reflection. This yields multiple 
scales in the convergence to the steady solution, and is one of the main causes of 
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the analytical difficulty of this paper. Eventually, we have the convergence rate 
(1 + + (1 — Q!)' 2 " of the boundary flux, Theorem [531 (1 —comes from the 

coefficient (1 — o) of specular reflection although specular reflection condition itself 
has no equilibrating effect. (1 + is from the diffuse reflection condition where 

the rate is essentially the same as before. However, we only have diffuse reflection 
for a multiple of a, and so the convergence to steady solution is slower than the 
case of the complete diffuse reflection. For instance, it starts to converge only when 
t > before that the solution is simply bounded. 

2.2. The global exitence for boundary flux. In this subsection, we prove the 
global existence for the boundary flux function j(y, t). It should be noticed that we 
may associate the boundary flux with the backward flow of particles. Once particles 
collide with the boundary, both diffuse reflection and specular reflection occur for 
the Maxwell-type boundary condition. Note that diffuse reflection is stochastic 
and see m for more details. In contrast to diffuse reflection, specular reflection 
is deterministic and has no equilibrating effect. In the following discussion we first 
give a solution formula of boundary flux for general domains. Here we assume 
temporarily that the accommodation coefficient is variable, 0 < a{y) < 1, for 
explaining what difficulties arise from this assumption. 

Fix y £ dD,t > 0, the boundary flux can be written as 


(32) 

j (-li •n(y))gi„(y-|it,4i)d|i 

+ J (-^1 ■n(y))a(y(i))MT(y(^,)(4i)j ^y(i),t- 


+ J (-^1 ■ n(y)) (l-a(y(i))) g ^y(i),t- 

l«ii 

= j®(y,t)+i^W(y,<)+F;W(y,t), 


where is a direct contribution of initial data, both and are the events 
that boundary collisions are more than once. More precisely, the first boundary 
collision takes place at y(i), D^^'> and E^^'> represent diffuse reflection and specular 
reflection of the backward flow respectively. We can continue to write down the 
formulas for and E^^\y,t): 


D^^Hy,t) = (y, i) + iy,t)+ -D^ (y, t ); 

E^^Hy,t) = (y, t) + (y, <) + E^ (y, t ), 

where 

Dil\y,t)= J (-Cl • n(y))a(y(i))MT(y(^^)(Ci) 

(-C2 ■ n(y(i))) 5 ™(y(i) - ^2)d^2d^i^ 
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(-€i • n(y)) a(y(i))MT(y(^))(^i) (-^2 ' n(y(i))) 


|y-y(i)l ^ |y(i)-y(i,i)l 
llil ^ IS2I 


^ ^ny . |y-y(i)l ly(i)-y(i.i) 

«(y(14))^T(y(,,„)(^2)j (^y(i,i)-^-- 


1^2 




DilUy,t)= J (-€1 • n(y))a(y(i))MT(y(i))(^i) (-I2 • n(y(i))) 

|y-y(i)k iy(i)-y(i.i)i 
' rein ^ IS 2 I 

, |y-y(i)l |y(i)-y(i,i)l 

(1 - «(y(i..,j s (^y(.,i|.< - —- 




Elu\y,t) = 


i -^1 ■ n(y)) (1 - a(y(i)) 




l«ll ^ liTl 


gin ( y(i) €1 (^ 


liil ) 




Edl}(.y^t) = 


(-^1 • n(y)) (1 - a(y(i)) 


|y-y(i)L |y(i)-y(2)l 
isii ^ m 


/ , |y-y(i)l |y(i)-y( 2 )l 

a(y(2))^T(y(,))(li)j (^y(2),^--— 




EiVeiy,t) = 


|y-y{i)l ^ |y{i)-y{2)l 

‘ -iiTi- 


1 — a 


(-€1 ■ n(y)) (1 - a(y(i)) 

y-y(i)l |y(i)-y(2)l 


9 y(2) 


t — 


111 


nil 


,1^ Mil- 


is the event that the backward flow reaches initial state after diffuse reflection 
occurs once. is the event that diffuse reflection occurs again at y(i 4 ) after the 
first diffuse reflection at y(i), and oipi is the event that specular reflection occurs 
at y(i.i) after the first diffuse reflection at y^^)- is the event that the backward 

flow reaches initial state after specular reflection occurs once. E^'j: is the event 
that diffuse reflection occurs at y(^ 2 ) after the first specular reflection at y(i), and 
Esll is the event that specular reflection occurs again at y( 2 ) after the first specular 
reflection at y^;^). Note that both and E^^ represent the contribution that the 

boundary collision takes place exactly once. More precisely, and E^^^ represent 
exactly one diffuse collision and exactly one specular collision, respectively. If we 
want to compute the contribution that the boundary collision takes place exactly 
twice, we need to take Eipi, E^^lj and E^ll into account. In other words, we 
must proceed to write down their formulas. Then there are four events arisen for 
exact two boundary collisions: {diffuse, diffuse), {diffuse, speeular), {speeular, dif¬ 
fuse) and {speeular, speeular). One can repeat this process inductively to compute 
the event that the boundary collision takes place exactly n times. In that case, 
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we need to handle 2" possibilities. That would make the solution formula lengthy 
and complicated. That is one of the main causes of the analytical difficulty of this 
paper. Another tricky problem is the variable accommodation coefficient a(y). In 
this paper we assume Q!(y) = a is a constant, this assumption not only makes the 
solution formula easier but also allows us to estimate all combinations of the events. 

We explain this by considering binomial expansion formally: 


n 

(33) {adiffuse + (1 — a)specular)'^ = {adiffuse)^ ((1 — a)specular)"' ^ 

k=0 


Then all combinations of events, R.H.S. of ([33]), can be dominated by 


{adiffuse + (1 — a)specular)" = 0(1)(q; + (1 — a))" = 0(1), 


for example, if we can show each term of diffuse and specular is bounded. In other 
words, we can treat the effects caused by diffuse reflection and specular reflection 
independently when a is constant. For the variable accommodation coefficient Q:(y), 
the problem is more delicate and might involve different techniques. This will be 
our another research work in the future. 

From now on we assume 0 < a < 1 is constant. We define the following notations 
inductively: 


y(o) = y(fci.....fc,.o) = y(fci,...,fco> 



11 


y(fei,...,fc,_i,i) - vb 

ii = 4r'-2(r' 




where y(ki,...,ki) indicates the location of particles via the backward flow process 
that: 


(fci — 1) specular —7> diffuse —>■ {k^ — 1) specular —7> diffuse —>■ 

• • • {ki-i — 1) specular —diffuse —>■ (fc; — 1) specular . 


According to the above discussion, we can find the solution formula of the boundary 
flux for general domains. 
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n { k k pi 

(34) j(y, ^) = II S H H J n ■ n(y(fci,....fe._i)) (1 - 


k —0 \ I —1 fci —I (^i^k) 

^(fc’i. 1=0 




\k-ki-...-ki 


9in y{ki,...,ki,k-ki-...-ki) ^l+l^ ‘ (j' HII 


i=l j=l 


|y(fci,...,fc.-i,j-i) y(fci,...,fci_i,j) I 


‘ |y(fei,...,fei,i-i) y(fci,...,fc,,i)l\ ^k-ki-...-k 


- E 


... d$,i 


(-^1 • n(y)) (1 - af 


|y(i)-y(i+i)l ^ |y(fc)-y(fc+i)l 
°<‘-£ —— 


(y,„ - &t - E 


n-\-l n-\-l « t 

7_1 7* I I 7, _7 ^ /_1 


Z—1 ki-\-...-\-ki—l 

^(kl,...,kl) 


.(-^.+1 • n(y(..,....M))) (1 - 

y{ki,...,ki,n+l-ki-...-ki) = k / , / , \ei~^\ 

i=l j=l I 

1 ) y(fei,...,fei,i) I (^n+1 —fci —... —\ 7jf 

TTTT, =^i+i «4;+i---a4i 

i=l l4;+il / 

[ ■ n(y)) {1-ar+^-g fy(„+i),t - " y(-+i) < ,g»+A 

, |y(i)-y(i+i)l ^ i=0 1 / 


i>E 

•i = 0 


where 


^(ku-.-M) ~ 

^ Y^^ |y(fei,...,fci_ij-1) ~ y(fci,...,fci_i,i)l ^ ‘ |y(fci,...,fc,,i-i) ~ y(fei,...,fc,,i) 

|eJ-i| Z^ |ei-i| 

2=1 j — 1 I 2=1 I^Z-1-1 I 


|y(fci,...,fci) y(ki,...,ki,i) 


if fc — fci — ... — fc/ = 0, 


|y(fci,... ,ki ^k—ki —.. . — ki) y {ki^...,ki,ki,k—ki — ... — ki + l)\ 


a k — ki — ... — ki >0 
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B 


{l,n+l) 

I ki 




i=l 3 = 1 

Note that 

•,ki + ...+ki—n+l 


1^^ 


i-11 


i=l 




(...) 


^k22.%'M + ---+kl=n+l 


I 

n • n(y(fei....,fc,_i)) (1 - . 

i=l 

( ^ ^^ |y(fci,...,fei_i,j-i) ~ y(fci,...,fci_i,j)l I 

J y(ku...M)^^-2^2^ -- I 


i=l 3 = 1 


le 


In the present paper we assume spherical symmetric domains and therefore we 
can make use of this symmetric property to obtain more precise formulas for the 
boundary flux by using change of variables: 


if d= 1, 


Si = 


|y(fci.....fci-i ,0) ~ y(fei,...,fci_i,i)l 


I [cos()ii = 


if d = 2. 


We define 


(35) 


For d= 2, 


ff(a) = (— ) e if d = 1, 

\ ( 7 / 

. 712 \ CT / 


n { k k pi 

d(y,0 = ^J^ ^ [ n(l-a) 

k—0 \ 1—1 ki-\-...-\-ki—l 


ki-1 


-AG.fc) 


i=l 


\k—ki — ... — ki 


, ./2i?T(y(fc ))si j2RT{y^^ 

-) - --(-«G. ^ n(y(G.*„)) 

I k—ki — ... — ki I I 

( ^k—ki — ... — kifi \ '' \ '' \'y{ki,...,ki,i—l) 1A ^k—ki — ... — ki\ 


i=l i=l 

d^ij^idsid4>i... dsid4>i 
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n,t |y{i)-y(i+i)l ^ ly(k)-y(k+i)\ 
- fTTI -<- TTS - 


ISII l«JI 

(-^1 • n(y)) (1 - fy(,) - 


1=0 


irii 


n+l n+1 


\ki — l 


aG 


EE / 11(1-«)' 

i — 1 fci+,,,+fci— i ^ 

®(fc’i,...,fc,) 

y2i?T(y(fc^_ y(2E7XyEEE 


ki 


L 

5(y(fci,....fe,)’i-E 


n+l —fci —... — ki 


, |y(fci,...,fci,i-i) y(fci,...,fci,i)I 

^ llUl 


i=l 

d^ij^idsid4>i... dsid4>i 


(-li+i •n(y(fc^ (l-a)"+i ' 

) 


.-ki 


n+l —fci —... — ki 

Z+1 


+ 


(-€i • n(y)) (!-«)”+ 5 y(„-Hi),< - 


|y(i)-y(i+i) 


l«il 


For d = 1, 


^n+1 

i=0 


111 I 


,^r 


j(y,o = ElE E 


(fel ,...,ki) 


Hd - .. 

. T ki 

t=l 

Z k—ki — ... — ki I I 

iy(/ti,...,fej,z—1) y(fci,...,fej,z)i 


-(yci.E 




)’ ^/-i-i 


k—ki — ...— 




d^i_^_^dsi... dsi 
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n,. |y(i)-y(i+i)l ^\y(k)-y(k+i)\ 

u<t- - 


ISII ' l«?l 

(-li • n(y)) (1 - ^ 1 d^i 


2 = 0 


111 I 


n+1 n+1 « 

EE / 

1 — 1 fcl + --- + fei—I (;_„ + !) 


n(l _ a)>-ag( .* !ll 


ki h 


L 


I n-\-l — ki—... — ki 

- 22 s,- 

i=l i=l 


Y- ly(fci,... ,ki,i—l) y(fci,...,fci,-t) I ^n+l — ki—. 

2. I^j-i| 


.-fc. 


...dsi 


+ j (-€i • n(y)) (l-a)”+^5 fy(„-Li),^ 

*>E 

i=0 1^1 I 


E 


2=0 


|y(i) ~ y(i-n)l 


j^n+1 

5 SI 


dll. 


To simplify the equations, we define 


(36) j-7(y>0 = 


(-^1 •n(y)) 


- 


l«f I 


t; I /-k/j. |y(d y(*+i)li ^k\ 

9in I y(fc) I I ^^17 


2=0 


1111 


(37) 

i?("+l)(y,t) = 


(-li •n(y))5 ^y(„+i),t-^ ^ ^ 


*^i5) 1^11 




(k) 

where jU represents the event that the backward flow reaches initial state after 
k times specular reflection and represents the event that the boundary 

collisions are more than n times and the first n+1 ones are precisely specular 
reflections. It should be notice that £^(”+1) is not the end. £(”+i) itself involves an 
infinite series and we will use the coefficient (1 — a)”’*'^ of £(”+i) to get the decay 
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for refined estimate later. Then we have for d = 2, 


n ^ k k p 

^ /(I-a) 

k—0 W=1 ki-\-k2-\ - \-ki—l 


k—k\ — ... — ki -(^“^1-^i) 

J in 


(y(ki,---,k,):i ~ ■si — ■ ■ ■ — ■si) 


Y[il - a)’^^-^aG 

i=l 

n+1 n+1 




dsd(t>t + (1 - o:)^jin\y^t) 


it-r± /fc-T-L p 

+E E - 

1 _1 I 7-. I I 7« _7 ^ 

y/2i?r(y(fc^_..._fe.))sA y/2^7X5EEE) 


I — 1 All “1-^2 ~l“ ■■■ ^ 


n(l-a)'=^-'aG I ct>, 

i=l 

and for d = 1, 


ki 


h 


-dsid(^i+(l-a)"+^£l(”+^^ (y, t), 


n [' A; A; „ 

•?'(y’^) = EiE E /(!-«) 

A;=0 I Z=1 A;i+A: 2-|-l-Ai^=Z 


k-ki-...-ki Ak-ki - ki) 

Jin 


(y(ki,---,fci)j ^ “ Si — • • • — S;) 


]J(1 - 

i=l 

n+1 n+1 


y2i?r(y(fe^... fc^))sA 


ki 


ki 


ds,+ (l-a)+|^^(y,f) 


J- / -L 

7_1 U I 7., - I \ U _7 

\/2^E>EEE)®A \J‘^RTiy(k^,- ,ki)) 


l—l ki-\-k2-\ -hA;i=Z' 




-ds,+(l-a)’^++(”+i)(y,t). 


Now we are ready to prove the global existence of the boundary flux function. 


Proof, of Theorem \2.2\ 

To compute the boundary flux j(y,<), we need to take all events into account. 
In other words, we have to sum up all events for each boundary collision. Hence, 
we have the following infinite series due to the above discussion. 


(38) i(y,t) = XI+1 (y^^)+ E E 


a 


,k—ki — ... — ki 


A;=0 k 

X J j 

0<si + .. .+sj <t 


Z—1 All + Ai2 + ■ ■ ■+A;^—Z 


(k-ki - ki) 


iy{ki,- ,ki)’^ — Si — ■ ■ ■ — Si) 


n (1 - a)'=-+G {^^i, d = 2 


0(1 -a)'=--+id 

i=l 


-dsi 


d= 1 


The index k here means the exact number of boundary collision and we will show 
the convergence of the series. From (1551) and noting that |^*| = |^| for each i, we 
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have 


^(1 - 0')’^jin\y^t) 

k^O 

oo „ 

(39) 


i=0 2=0 


- $. ■ n(y) 




(1 + C)^ 


drjd^ 


^ IlSirallc 


-l-nCy)! 


J (i + C)'' 

= ^^(1) ll5in|loo./i 


dc 


With (l39ll . we rewrite (l38ll as 


OO ( k 


j{y,t) = 0(1) ||g„|| 


oo,/j, 


E E E 

k—0 I 1—1 ki-\-k2-\ - \-ki—l 

- 


.(fc-fci- ki) 

-Jin 


0<si-\-... + si<t 


n (1 - a)'=^-iaG ^ = 2 , 

0(1 -Q)fe^-iag 


00 00 


(40) =0(1)||9.„|1 


St—i; (1 _ 


00 ,/X 

Z=1 fei = l = l fe=/ci + ...+fei 


.(k—ki - 

Jin 


0<si + ...+Si <t 


0(1 - a)^^H - Vg^ y^M, ■>=.)) 


By plugging (IMl) into (HHl) and direct computations, we have 


J(y,t) = 0(1) ||g„|l^ 


00 00 00 

i+Di^tE-E 

. 1—1 ki — 1 ki — L 


O'^si-j-.. .-j-si <^t 


0(1 - a)^’G d=2, 

0 (1 - a)^'H - 


(fcl,--- ,ki)) 


dsi ^ d = 1 







































20 


HUNG-WEN KUO 


= 0(1) ||5i„|l 


00 ,/J. 




'1 — a' ^ ki 

1=1 fci=i ^ 


ki = l 


h 


n 


(^||G||^„y2i^)ds„d = 2, 


0<Csi-|-...-)-Si'<t 


t\\ {\\H\\^^V2m^)ds. ,d 


= 1 , 


00 


I 1^1 \fcl=l 


{7v\\G\\^^V2m^ty 

V _ , a — z 

{WHw^^vmf^ty ^ 

,, , a — i 


It is easy to check that 


and therefore we have 


^ (1 - a)'^^ 

-r-= -Ina, 


ki — 1 


ki 


j{y,t) = 0(1) ||5„|l 


I e^r^^l|G|hooV2lj7^i^ ^ ^ 2 , 


□ 


2.3. Preliminary Estimates. The discussion of this subsection applies to all 
space dimension. To avoid complication in notations, we treat only the 2d case. 
And we will use the following Law of Large Numbers to get a refined estimate. It 
can be proved by the similar argument as in m- Therefore we omit it. 

Theorem 2.5 (Law of Large Numbers). There exists some constant G > 0 such 
that, for any 7 and m with 7 /(mn)^+T > G, 

log (7 + 1 ) 


J Hn{a)da = P {^ < l^i+- • .+A„-nE(Xi)|} = 0(1)- 

,X<|cr-nE(Xi)| 

where Ai, A 2 ,..., are i.i.d. with the probability density function 


yd+l 


H{a) = 


*/d=l, 

J G{(j),a)d(j) ifd = 2, 


and Hn is the probability density function of the sum of i.i.d. random variables, 

Ai+A 2 + - • ■+Xn. Notice that iL„ is a convolution of H itself, Hn = {H * ■ ■ ■ * H). 

^ ^ 

n times 

In order to get a refined estimate, we start with the specular reflection: 

m— 1 


(41) j{y,t) = 51 “ oifj\n{y^t) 




t/2 t 


k^l 


+ a55(l-«)'' ^ J J j{y(k),t-s)G{(f, 

-7r/2 0 


^\J{yf^k)) Y^2i?T(y(^.)) 


dsd(j) 


+ (l-a)™E(™)(y,t) whend = 2 . 
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A;=0 


1 _ 1 ^ 


k^l 


k ' k 

+ {l-a)’^E^"^\±l,t) whend=l. 


Define 


m— 1 

A;=0 


and rewrite dm) as: 


(42) j{y,t) = + (1 - a)-i?W(y,i) 


7I-/2 t 


+ 0^(1 


/c=l 


,t- s)G'((/), 


-7r/2 0 


dsd(j) 


By iterating n times, we have 

(43) j(y,i) = 4”^^(y,t) + (1 - a)'”i?(’”)(y,i) 

n —1 m m 


2=1 ki — 1 


ki^l 


0<CSi+S2H- 


SiY^2i?r(y(fe^^... fc.)) v^2i?T(y(fc^_... ^,^)) 

11 G(0j, 7^- ) 


i=i 


% 


{•^ir^(y(/ci.-.fci)>i- Si -Si) + (1 - a)™£^^™Hy(fei,...,fe,),i- Si - Si)^ ds,d(l), ■ 

mm « 

+a” E/ j(y(fei,fe 2 ,-.fcn)>^-si-S 2 -S») 

7, _1 7^ _1 


ki — 1 


fen = l 


0CSi+S2H- 


l[G{cl>i 


i\J‘^RT{y{ki,-,ki)) Y^^2i?TljEEfci))' 


dS’fidcjyYi * * ■ ds idc^i. 


2 = 1 


To estimate (Pi) . we define 


j‘lk^l.^k,)iy^t)= J 

0<Csi+S2H- ^ ^ 


* Sjy'2i?r(y(fe^_... y/2i?r(y(fe^_... 

11 ^(Qj! )■ 


{'^ir^(y(fei,... i - Si - Si) + (1 - a)™-E(™)(y(fe^_... j,.), t - si -s,)} ds,d((). 


• • dsid^i 


• • (isici(()i 
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and 

'^(r;^.fc„)(y>0 = J Si- S2 -s„) 

0<Si + S2H- 

" Si G^2i?T(y(^.^ ... 

J| G((;!)j,- - -- - - dsndcjin ■ ■ ■ dsid(j)i. 


Recall 


m— 1 


•4n (y> (y> t) 


k—Q 

m — 1 

= 0(l)^(l-af ||5»IU 

= 0(1)to||5„|1^ 


Moreover, we have for i > 1, 

jln\y^t)= [ 

-<ICI<- 

= 0(l)llff^n|lc 


-I • n(y)5*„(y(fc) -^'"(t- kh), 


fc|y-y(i)l ^ 1,1 ^(fc-M)|y-y(i)l 
t ■^lsl< t 

d+1 


= 0(1) ||5i„|l 


((fc +i)|y-y(i)l) -(fc|y-y(i)l) 


. d+1 


oo,^ 


^d+1 


= 0(1) ||5i„|l 


2k+1 

~1^ 


for d = 1, 


1 3k'^ + 3k+l 
ts 


for d = 2. 


Thus we have 


m— 1 


(44) 


(y> t) = 0(1) lldmlloo,^ H 


fc =0 


2/c-l-l 
3fc^ + 3fc + 1 


i3 


for d = 1, 

for d = 2, 


= 0(1)||5™||, 


1 -I- (1 — a)m‘*+i 


In order to estimate the remainder terms, we define a priori bound of boundary 
flux j{y,t). 


for d = 2, 


Definition 2.6. Define the a priori bound J" by 

_ I (ii^ii-r) (*> 

J (t) = sup < 

[ ( b+(s)| + |j-(s)| j for d = 1, 

where j±(t) = j(±l,t). 

From (Ei) and (EH), it is easy to show 

F;(-)(y,t) = 0(l)(:7(t) + ||5.„IL,J. 
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For i < n, we divide into the following two parts: 

‘iL,*L) ‘~~ 






1=1 


_.(i,m)sloWf ,x , ■{z,mjrare, ,s 

where 

fclCTl 


• (z,m)r 


= < 0 < 


... + 


< 


'T, t 


■ ■ ■ ^2RT^kk,...M) V T* 2 j ’ 


^2 = < < 


t k\(T\ 


kiCi 


T*2 ^2RT^kk,...M) 


<t\. 


For (y, t) , the time needed to trace back to an interior point is at least 


t/2: 


ki<7i 


+ ...+ 


kiCTi 


< 


'T* t 


^2RT^kk,...M) V T* 2 

^ _ / kiai , , fcjCTi 


I y^2i?T(fej) y^2i?T(fei 


> t - 


t ^ t 
T*2 ~ 2' 


Thus 


.{i,m)slow, N 


< / sup 




'y<S<t 


J, 


(m) 


+ (l-ar 




(s) X / JjG((/)i,cr;)d(/)i(icr/ 




(45) 


Note that the estimate (H51) relies merely on the smallness of the speed. 

For t), the time consumed to trace back is at least ^J^t/2. There- 




fore, 


' T* t 
T* 2m 


\/ 2RTf < (Ti + ... + CTi, 


.(i.rn \‘rare. / , \ 

<0(1) (m ||5*n|loo./x + (1 - a)"" (llff™|loo,^ + ^(t))) 




2^ \/2/tT^+ C^O"! -j-.. .-1-fTi 

=0(1) + J^(t)) X 


Hi{a)da 


(46) =0(1) (m + J{t)) X Pr |xi + ... + X, > . 
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Note that n, m, the index of ^.) (y) i); are variables at our disposal. Through¬ 

out this paper we assume t/mn ^ 1. Recall that our final choice of n and m is 
mn = r S (0, (d + 1)“^)- Thus r < 1 and so, for large t, mn ^ t. Since 

E(Xi + ... + Xi) = iE(Xi) ^ i < n (H51) represents the prob¬ 

ability of a rare event. We now apply the law of large numbers, Theorem 12.51 to 
estimate (PI): choose the truncation variable 7 to be V2i?T,|. Under the assump¬ 
tion t/mi > t/mn ^ 1 , we have 


IT, 


t 


(mn) 'i+i (mn) < 1 +^ 


— > 1 . 


c {7 < k-^E(Xl)|}, 

Therefore, we can apply Theorem 12.51 to obtain 

t A' 7 \ .{i,m)raref r^t^ \ / II II i ™ ^ * log(t-(-l) 

(47) = 0(1) [m\\g,n\\^^^ + J{t)^ --. 

Note that the estimate (H71) relies merely on the law of large numbers. 
For ), we conduct a similar decomposition: 

4i’ ^ - ... - 

n 

^ Y[G{(l)i,ai)d4>idai 

,m) / ,\ , T('n,m)rare / ,\ 

h\(Ti kjiCji 


i=i 


where 


= <^ 0 < 


-\/ 2 RT(^ 


< 


r* t 


52 = 


r» t 

< 


fclCTi 


k2Rr(fe^.....fc„) V T*2 

klCTji 


...+ 


<n. 


T*2 y'2RT^k„...M 

-TXT ^ , r ■(n,m)rare , r(n,m)rare , i . • 

We can apply the same argument ot to 7)^.^ to obtain: 

fAO\ -i(n,m)rare, ^ lo§(l T 1) , ,/ ^ i 

(48) J(ki,-,ki) = 0{l)J{t) --, whenever t/mn » 1. 


We omit the details. 

Lemma 2.7. 

m—1 m—1 

^ ... 

ki—0 ki—0 

n—1 m—1 m—1 


a 

n—1 


E E ■ ■ • E (1 - ^ (1 - (1 - arr = o^y 


2=1 fei =0 ki—0 


2=1 


Since 


Jiy,t) =J<f’(y,t) + (l-a)™U™)(y,t) 

n— 1 m m 

+E“' E (1E(i-«)'''”'4fcE.fc.)(y’^) 


2=1 fci = l 




+a" ^ (1 - • • • E (1 - 

/;:i=l Ain —1 

putting (HSl) . (HTl) . and (H51) together we have: 
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Theorem 2.8. For t/mn ^ 1, 

=0(1) + (1 - o)-) 

m m 

+a" ^ (1 - • • • E (1 - 


fci =1 




where 


(49a) = J j - si - ... - s„) 


0<Sl+ .. . + SjT, <c Y ^ 


X . 


when d = 2, 


(49b) A|”E ,,^)(±l,t) = J j Si - ... - s„) 

0<Si + ... + Sti<^ ■§■ 


n« 


ii,^2RT(±l(fc^.V2j^r(±l(fc^^...^fc,)) 

h 


fc,-y - ^ - dsi J when d= 1, 

• (n,m) 


With the aid of the law of large numbers, we have estimated and 

remaining term ^ ^ consists of the main event, which requires 

more effort to estimate. 

To show the convergence of boundary flux j, we need to use the crucial conser¬ 
vation of molecular number, dMl) and (E51) . as 


dKd^ 


j{y.t) = ^ ini / 5 (x,|,t) 

'-'S\F>\ JdxR'J 

[ (j{y,t)s{x,^)--^g{K,^,t))dxd^ 

JDxRrf ^ US / 


1 

^1 

1 


(50) 


\D\ 

“'I€l< 

l>=-X(l)l V 

t 

1 

f 


1^1. 

/|=c-x 


1 

f 


1^1. 


l*(l)-i'(2)l 

’ l€l>- 

log(t-fl) 


Cs' 


log(t+l) 


(j(y,t)s(x,|) - ^ 5 (x,|,t))dxd| 


(j(y,<)s(x,|) - ^ 5 (x,|,t))dxd| 


= jm(y,t) +jmid(y,t) +jfiiy,t). 

It is easy to see that if we choose K = K{x,^,t) such that 

t |X - X(i ' 


a:- 1 < 


<K, 


log(t-bl) |X(i)-X(2)| 

then we have 
(51) 

|x-X(i)| -t (AT- 1 )|X( 1 ) -X(2)| |X(1)-X(2)| |x-X(i)| - 1 -A:|X(i) -X(2)| 


t 


< 


log{t + 1) 


< 


t 
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Since the domain is spherically symmetric, it is easy to show that 

|X(1)-X(2)|“ 

and therefore K Ri iog(t+i) ■ With (|5T]) . (1^ . and (1501) . we have 


(52) 


Cs\D\ 


'l«l< 




oo 

^ i —1 


i -1 




27T 


((Rr(x(i)) 








(53) 


K 


jmid{y■> ^ 


27T 






2=1 


X MT{jc^k+i))ii) - 9 in{X(k) - l*(i -h- ... - tfe),|^)^ I 


2 tt y 

i?T(x(fc+i))y 

dxd^, 


(54) 


jfiiyd) = 


27r 


Cs\D\ 


[ ^ [j{y, t) - - h - ...-ti)] 

1^1 > log(t+l) ( i=l 

27r W 


j + <' - (“g<‘ - « (Tmwi 


^'r(x(K-^i))(€) -5 (x(k) - - ... - tK),l“) 1 Wxd^ 


K\ 


where 


. , |x - X(1)| + {k- 1)|X(1) - X(2)| ^ 1^, |x - X(1)| -h A:|X(1) - X(2)| 

Ak = s -7-< l€l <-:- 


Each component of j can be estimated in terms of J{t) and the fluctuation of 
j. Therefore, it suffices to consider the fluctuation of j. For t' < t, mn/t' <C 1, 
Theorem 12.81 yields 


00 ,/i 


(55) i(y,() = 0 ( 1 ) ( '"("‘"*''),y' + l) +(l - a)”‘) Ilfcll 

+o(i)( (’"")'y^i) +(i-°r)j(*) 

m m 

+a” i: (1 - »)*■-■ ■ ■ ^ i: <1 - “)‘"^‘(''S;:"h.,(y.') - A<r;:"d,,(y',t')). 


ki — 1 
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And note that a" Sfc^=i(l ~ ^ ‘ ~ ^ — 1- Therefore, we need 

only to estimate the fluctuation of ^ ^ and show that they are uniform for 

each (fci,..., kn)- Since 


(56) A!E:‘h.)(y.‘)-''S::h.)(y'.*') 


we may consider temporal and spacial fluctuation separately. We study the tem¬ 
poral fluctuation in Section 12.41 and the spacial fluctuation in Section 12.51 


Remark 2.9. For the Maxwell-type boundary condition, which is a convex combi¬ 
nation of the specular reflection condition and the diffuse reflection condition, the 
intricate dependence on the accommodation coefficient a yields serious analytical 
difficulties beyond those in [1^ and HU. One needs to consider all the events of 
particle colliding with the boundary many times, psp . In that case we need to deal 
with the problem that the diffuse reflections are coupled with specular reflections. 
Roughly speaking, for the events that specular reflections are more than diffuse re¬ 
flections, we need a new idea to obtain the decay rate even if the specular reflection 
itself doesn’t have the equilibrating effect. We achieve the aim through (l4^ and 
(HU). For the events that diffusion reflections are more than specular reflections, 
we may modify the analysis from the previous works m and m to get the decay 
rate, (HSl) . However, the appearance of specular reflection will slow down the decay 
rate. Finally, we succeed in combining all events via Theorem 12.81 


2.4. Temporal Fluctuation estimate. In this subsection we consider the tem¬ 
poral fluctuation. Recall 


. (n,m) 


{y,t) = 


/ 


0 '<Si+S 2 H-l“Sn<^ 



jiy{kxM,- ,fen)’ ^ — Si — S2 — • • ■ — S„) 


2=1 




)) 


ki 


ki 


~dsjidcj}^ * * * ds\d(^\ 


o< 


fclo-1 , 


fciCTi 


<t/2 


3 I y(A;i ,^2 ,••• ?^n) ’ 


t — 


kn (Jn 


n 

G{(j)i, ai)dand(j)n ■ ■ ■ daid4>i when d = 2, 


i=l 


^(n,m) 


(fei,... ,fc„) J j ,kn}y^ '®2 ■■■ Sre) 

0<Sl-t-S2-|- \ j 




i=l 


ki 


ki 















28 


HUNG-WEN KUO 


0< 9PT --1-I- 9BT -<*/2 

2R^(fci,fc2) ' 




klCTl 

2-RT'(fci) 


!(fci,...,fc„) 


H{ai)dan ■ ■ ■ dai when d = 1. 


We note that the kernel H{a) and G{(j),a) are smooth in u, and hence we may 
differentiate ^ ^ with respect to t directly to obtain an explicit expression. 

Lemma 2.10. Let n be any positive integer. ^ 'jiUjt) is with respeet to 

t. Their derivatives has two parts: 


{ki ....,kn 




The first term ^ is the boundary term: 




siy/2RT^k.^ . k^) \ yJ‘^RT{ki,...,ki) 


Sl+-.- + S-n-J :p-| 


X J I y{k^,...,k^)fi- Y 7^2 I when d = 2, 


(fcl,...,fcr 


,(±l,t) = - 1- 


siy/2RT(^ki,...,ki) 

ki 


/2RT(^ki,...,ki) 

h 


Sl+... + Sn — A 


X j I - Y —- I dsi • -dsn-i when d=l. 


The seeond term ^ ^ is the volume term: 


(v T-- 


* J 

Si + ... + Sr,=t-S 


ki ] ' 


^y2RT(^ki,...,ki) dG I ^ \ r< ( M s„y^2i?r(fcj ,,, 

X J ^ I J I X • • • X G I J 

ki u(j \ ki } \ kfi 


X dsi • • • dsn-i j{y(k^,...,k^),s)ds- 


f {ki,...,kn 


-d^4> when d = 2, 
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(60) = 


I X 




h 




ki 


da 


h 


X dsi ■ ■ ■ dSn-l s)rf5 ^ {ki,...,k„) d = 1. 

Ki kji 

The lemma can be proved by a similar argument as in mm and we omit it. 
The boundary term Bn can be easily bounded as following. First, we have 


’ ’ i<s<t ^ ^ ^ 


no (■>., .‘■d . i *, ■ ■ ■ 


Si+.-. + Sti — w 




A:/ 


0(1) sup (ibllicx,) (s) 

4<s<t ^ ^ ^ 


'(fci) 


-I-...-I-g 


JJ G ((/>;, a/) dsi ••• when d = 2, 


z=i 




l^(fe(™hn)(=^l’^)l = *^(1) ^'^P (b+(«)l + l•)-(■S)l 

^<s<t ^ 


Sl + --- + Sn — 4 




h 


ds\ • • • dSn—l 


0(1) sup (|j+(s)| + |j_(s)| 


k<s<t 


X 


H (ai) dsi ■ ■ ■ dsn-i when d = 1. 


: + ■■■+5 


1=1 


J'(fcl) ’ ■■■ ’ 2HT(fc^. kn) 2 

Hence given t' <t with t'/mn ^ 1, by the law of large numbers, Theorem 12.51 

'■* ' " m?rB' log(t' + 1) 

d3 ’ 

when d = 2, 


(62a) ,^^)(y,s)|ds = 0(1) ^sup (IIjIIl-) (s) x 


(62b) 


/ = ^(^) ®^P (l■l+(■s)l + b'-(s)l) ^ 

dt' ’ ’ T:<S<t ^ ' 


m^nlog(t' + 1) 

d2 ^ 


when d = 1. 
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We next turn to the major term Vn, First, as above, we have 
1 


(63a) )(y,t)| < - sup (||j||ioo) (s) 




1^1 


kl KT^7-±/ 

(63b) )(±l,t)| < i sup ( \j+{s)\ + |j_(s)| ) 


(PacTcj) when d = 2, 


i<s<t 


. 


1=1 


kl da 

To estimate dMl), we use the following lemma. 
Lemma 2.11. For any integer n > 1, 


dFa when d = \. 


(64a) 




1=1 


h 


dJ^add^ft) 

= o({nlogn)^ ) , 


(64b) 


i: ydFEFFhiHM- ■ ■ HM 


1=1 


kl 


da 


d^a 


= O ((n)2 ) , 


Consequently, 


j<S<t 


=0(1) ^sup (IIjIIloo) (s) X , whend=2 


=0(1) ^sup (|j+(s)| + |j_(s)|) X , whend=l. 


i<s<t 


Proof. 

Note that fc; > 1 for each I, so 


f: . ^ ^ 


1=1 


kl 


< 


da 
dG. 


_^ ! - uij 

^ ^ Y'2-RT'i^/j.^ J G(01 j CTl) ■ ■ * ' ' ‘ (7n^ 

1^1 

Then we follow the Lemma 3 in m to conclude the proof. 

The following theorem follows from Lemma [2.111 together with (IS^ . 
Theorem 2.12. Let t' < t, then, for 1 <C t'/mn, 

Gn.m) , . _ . (n,m) / ,/x 


□ 


= 0(1) ^sup (llillioo) (s)x 


log(t'+ 1) /log? 


■V <s<t 


t'3 


(t — /) ) , when d = 2, 
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A (n,m) 


(±l,t) - A 


(n,m) 

(ki,...,kn) 




= 0(1) sup (|j+(s)| + |j_(s)|)x 


! w?n\og{t' + 1) 



when d = 1. 


From Theorem 12. 121 (1551) and (1551) . we obtain: 

Corollary 2.13 (Temporal Fluctuation Estimate). Lett' < t, then, fori -C t'fmn, 


t') 


0(1) + (1 _ o,™) (||,.„||^^ + J,„) 

1 

+ 0(1) sup (|IjIIl“) (s) X (t-t'), for d=2, 


j{±l,t)-ji±l,t') = 0{l) 


m{mn)^ log{t' + 1) 
+2 


+ il-a)A (||5™IL,^ + ^(t)) 


+ 0(1) sup (|j+(s)| + |j-(s)|) X (t-t'), ford 


'= 1 . 


2.5. Spacial Fluctuation Estimate. In this subsection, we investigate the spa- 
cial fluctuation. 

Theorem 2.14 (Spacial Fluctuation Estimate). Suppose that t/mn, t/mN, 

1, 0 < g < 1, y, t/ S dD. Then 

= 0(1) ( (("")'+ + (1 - a)") (b.„ll„^, + J(1)) 


^<S<t 


+ 0(l)^sup (IIjIIl-) (s) X I f mN + 


logN 

N 


when d = 2, 


ii(+i,i) - = 0(1) + (1 -o)"*^ (j(t) + iis,.ii„ J 

+ ((;) ’ '• + ( r )') fz., (^ 

when d = 1. 

We consider the one dimensional case first, which is much simpler than the 
multidimensional cases. 

One DIMENSIONAL CASE, d = 1. 

This case differs from the multidimensional cases in a fundamental sense: unlike 
the multidimensional cases, the boundary, comprising of two points, is discrete. 
This makes the one dimensional case much easier. Instead of processing we 
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directly estimate: j{+l,t) — Recall 

m— 1 

j(+l,t)= ^ (1 — 




0 


fc=l 


From 


f ^ s^2RTii-m ^2RT{i-m ^ 

J k k ' 


we have 


.,-M) = 




Now consider 


m—l 


j( + l)0 - - Oif j\n { + 

+of (!-„)- / (;„-I)d* - .) - ;,-M)) 


fe=l 


— a 


{ k k 

+ (l-ar {E<^'^\+l,t)-j{-l,t)). 

As before, 

m — l 9 


(t + l)2 


fc =0 


(1 - a)™ [E^-^){+l,t)-j{-l,t)) = 0(1)(1 - a)™ (j{t) + ||5i„IL J 
and we decompose 




k ' k 

ft t/2 


0 




t/2 
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where 0 < <7 < 1 is to be determined. And then 




k—1 i 


ds 


= omnhY^ii - 




t 2 


m « 1.2 2 

■ ^(1 - af-^ /(•••)= 0{l)J{i)aY,{l - af-^- = 

k—1 ^ k—1 


t/2 

i/2 


2 




^(l-a)'= (^...^ = 0 ( 1 ) ^sup (|j+(s)| + |j_(s)|)a^(l-a)'= 

k — 1 2 

sup (ij+(s)i+ij_(s)iy 


m'' 

^<s<t 

Finally, we consider 
to 


_i)>, t _,) _ ,(_i,,)) H(»w«n(-i)‘), 


0 




Lo H, Vises, AIM 


and it is easy to see that 






ds 


< sup |j'(±l,t) - j'(±l,t')l 




■Bi-A- / (i((-i)L,)- j(-,0)AIM,. 




= a 




fc=l,fc even 


< X—^ (j'( + l.i) - j(-l,i)) ■ 

Z — Of 


Hence 

1 


2-0 


(ji+l,t)-j(-l,t)) = Oil) ^— + (l-a) 

+ sup \j(±l,t) - j(±l,t')\ 


hv 




ds, 
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and it follows 

|j(+l,t)-j(-l,t)| = 0(1) +1) + ^ + (1 - a)™) [j{t) + ||5„IL,^ 

+ 0(1) (^)' sup (b+(s)| + |j_(s)|) 

^ ^<s<t ^ ' 

+ 0(1) ^sup ( |j+(s)| + |j_(s)| ) X 

Subsequently, we can apply the temporal fluctuation estimate, Corollarv l2.1111 to 
the term sup^_^g^^,^^ |j(±l,t) — j(±l,t')|. Hence Theorem [2+4] for d = 1 follows. 
Note that t — > | so long as t 1. 

Multidimensional cases, d = 2. 


As noted before, to estimate spacial fluctuation we invoke another variable N. 
From now on N will be the index of . 

The boundary dD is unit circle, so we parametrize it by the polar coordinates. 
Given two boundary points y and y', let y' be point of degree zero, and denote the 
polar angle of y by 9. Denote the relative polar angle of with respect 

to Then for 1 < j < ki, the relative polar angle of y[kx,...,ki-i,i) 

with respect to y[ki,...,ki-i,i-i) is also 9i because of the specular reflection, i.e. 9 + 
ki9i +.. . + ki9i stands for the absolute polar angle of y(ki,...,ki)' ^i^ace dD is the unit 
circle, 9i = tt— 2(^;, To simplify the notation, put T^ki,...,ki) = T(9+ki9i +.. . + ki9i), 
and = T{ki9i + ... + ki9i). Under this coordinate system, we have 


(65) 






N 






1^1 


xj(9 + ki9i + ... + t - - ... - - ] d^ad^4> 


— / - +----I— / 

\/2-R^(fci) y/^^'^(ki,...,k 




n) 


/ T* t 
2 


1^1 


X j 


+ ... + k]\f9]\f 


k-\ (7-[ 


kw<yN 


^2Rf^kg,...,k^) 


d^ad^<i), 


where T[ki,...,ki) ~ T(d + ki9i + • • • + kpfdN — ^)- 


Remark 2.15. In [TU] and m we can deal with the case of spherical domain in 
for d = 1,2,3, but in this paper we only consider the the spherical domain in 
for d = 1,2. As in our previous works m and m we need the symmetric 
property of domain to calculate exactly the spacial fluctuation. The symmetry of 
the boundary allows us to tract the exact location of the particle after multiple 
reflections. This is an essential ingredient of our analysis on treating spacial fluc¬ 
tuation. For two dimensional case, thanks to polar coordinate we are able to tract 
the exact location of the particle on a circle after mixed specular reflections and 
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diffuse reflections. This allows us to conduct a change of variables, (IS51) . and to 
estimate the spacial fluctuation, (Ei)- However, in three dimensional case there 
is no universal coordinate to tract the exact location of the particle on a sphere 
after mixed specular reflections and diffuse reflections. The three dimensional case 
might require mathematical analysis different from ours. We hope to return to this 
problem in the future. 


From (ESI), we have 
( 66 ) 


N 




=+...+ 


'(fcl) 






Ifc 

1=1 


9 


2kiN 


cri 


X j kiOi + ... + kNON, t - -... - , - d^ad^cj) 


N 






= + ...+ 


^ /2RT' " , /2RT' T T 

V (^l) V 


X J ( fci 01 + ... + - ... - -) d^ad^cji. 






2RT' 


The two terms in (I66|) differ from each other in three places: the domain of in¬ 
tegration, the angular variable of the transition PDF G, and the time variable of 

j- 

We now break the spatial fluctuation into three parts: 

A!f;':b„,(y.*)-Aif;:',.„,(y'.()=A|"':i,„,(9,o)-A<»;y,„,(o,t)= um + u ,. 

where 

/ \ 

Ui = 


■" ..Aa,) 


(ki,...,kj^) 








2RT( 


(ki ,...,kj^) 


N 


U2 = 




=+...+ 




^ /2RT' ” ^ /2RT', T T 

V (*=l) V (ki,...,k^) 


kN<yN 


j I ki9i + ... + k]\[9]\f , t -y . - , 


— j ( ki9i + ... + t --... — 


kN<7N 


4 '- 


(fci) 


J2RT' TT 


d^ad^cj), 


d^ad^cj), 
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N 


N 


Us = 






, /2RT' ! " 


T* t 
2 


IlG 

\i=i 


2kiN 






Z=1 


X j [ ki 6 i + ... + UnOn, t - - ...- 


kNf^N 


d^ad^cj). 




(fci) 


^2RT, 




Ui,U 2 ,U 3 register the difference in domain of integration, the time variable of j, 
and angular variable of the transition PDF G, respectively. 

We now proceed to estimate Ui, U 2 , and U 3 . Consider first Ui. As noted before, 
Ui registers the difference in domain of integration. Denote by A©i? the symmetric 
difference {A\B)iJ {B \ A). Since ki <m for each i and T > T* on the boundary, 
one can observe that both of the events 


S'x — 


<02 = 


kidx 


+ ... + 


kN<JN 


2RT, 


(fei) 


< 


2RT, 




kiai 




+ ... + 


kN<JN 


< 


contain 

So we have 

(67) 


cri 


+ ... + tr^v < \/ 2RT^ 




T* t 


T, t 
T *2 

I T, t 
T *2 


TL t 


A ©<^ 2 © •^cri + ... + (T 7 v © \/ 2 .RT^\l r ■ 


This implies that A © <^ 2 ) and thereby {7i, is a rare event, and can be estimated by 
the law of large numbers, Theorem l2.5l 


Lemma 2.16. Assume t/{Nm) ^ 1, 

( 68 ) 


u, ,0(1) "“"^ y + 1) , Jfty 


Proof. 

From (IHT)) . 


G{(j)i , CTi) • • • G{(I)n, (TN)d^ ad^(l). 


\Ui\<J{t)x2 J 
Applying law of large numbers Theorem 12.51 with 7 = y/2RTl ^^|, we conclude 

(IS51) . □ 


U 2 records the difference in time variable of j, which is exactly the temporal fluc¬ 
tuation of j. Therefore, we can apply our previous estimate in temporal fluctuation 
to this part. 


Lemma 2.17. Suppose that t/mn,t/'knN, N ^ 1. Then 

{mn)^ log(t -I- 1) 


U 2 = 0 ( 1 ) 


(J +(!-«)") (llfclU 

0 ( 1 ) s„p (|U||,.)(,)x''‘‘>'" . 


y<S<t 


N 



( 69 ) 
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Proof. 


N 


U2 = 




9 


2kiN 




1=1 


. 


kN<yN 


j ki9i + ... + ki\[ 0 ]sf, t - , — ... - . , 


fc/VO’AT 


+ + .. 


d^ad^c 


( 


*^1 _|_ j_ '-iv " i\ _^ 


\ 


|cri + ...-|-crN-ArE[Xi]|>A J|cri + ...+crN-AE[Xi]|<A 
k N^N ^ . / T, t 

kjw) 


(•••) = c /21 + U 22 , 


For U 21 , since 


fclCTi 


kN<JN 






T.^t 

< Y ^2 ^ 


fciCTi + ... + kfqaN < '/2RT*\ 
we have 


' r, i 

7^2 


•C< 


kiai 


kN<JN 


2RT, 


(ki) 


< 


2RT, 




IC/ 21 I < ^sup (||j||i^)(s) 

J G{(j)i,ai)---G{(j)N,crN)d^ad^(j). 


^<s<t 

X 2 

|o-i + ...+crN-IVE[Ai]|>A 


Applying Theorem 12.51 with 7 = A^, we obtain 

(70) C /21 = 0(1) ^sup (IIjIIl-) (s) X 


i<s<t 


logN 

N ■ 


Note that the prerequisite N/Ns ^ 1 of Theorem 12.51 is satisfied since iV >> 1. 
Next, for C/ 22 , by Corollary 12.131 

IO 22 I 

<sup||j(y,s} - j(y,s')| : s,s' e (^t - mN ^ , y e doX 

(71) , 0 , 1 , ^ ^ 

1 

+ 0(l)^sup (|I/IIl~) (s) X 


From dZOl) and (ED we conclude dUD- 


□ 
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Finally, we investigate C/ 3 . C /3 involves only the angular difference of those PDF 
G. No difference in boundary temperature are included. Therefore, the estimate of 
C /3 is reduced to the constant boundary temperature case, as in [10]. Since all the 
9 dependences appear only in the N copy of G, C /3 is a function of 9. Moreover, 
= 0. By direct computations. 


(72) 


dC /3 

d9 


< sup (||j||i=o) (s) 

|<s<t ^ ^ 




N 


1 dG 




1^1 




The RHS of (17^ can be derived by the similar argument as in m- Thanks to 
ki > 1 for each I, the following lemma allows to obtain a decay of C /3 in N. For a 
proof, see m- 

Lemma 2.18. 


N 


dG 


^ G((/)i, CTi) • • • ■ ■ ■ G{(j)N, ctn) 




d^cj) 


o((/VlogiV)5), 


Therefore, 

(73) I C/ 3 1 < 


p'K 

dUs 

L 

d9 


d9 = 0{l) sup (lllllioo) (s) X 

i<s<t “ 7 


logN 

N 


Under the assumption t/mn,t/mN, N 1, patching (1551) . (IMl) . and (1751) to¬ 
gether we have 


(mn)^ log(t + 1 ) 
t3 


{N,m) 


= 0 ( 1 ) 


+ (l-a)-] (lls.nlU,+ //(/)) 


+ ( 


0(1) ^sup (IIjIIl-) (s) X ( 

2<S<t y 

Plugging this to (|55]l . we conclude Theorem 12. 141 


/logN 


logn 


mN . 


2.6. Convergence of Boundary Flux. In this subsection, we prove our main 

theorem. Theorem 12.31 for free molecular flow. 

To apply a priori estimate, we need to establish the boundedness of j first: 

llj’llioo = 0(1) llsmlloo fj,’ cf. Lemma [2.211 We apply a priori estimate twice: in the 

first time we obtain a rougher estimate, the boundedness of j, and in the second time 

1 

we use the boundedness of j to obtain the convergence rate {at + 1 )“^^ + (1 — a)‘™ 
of j. 

Now we recall j(y,C) = jiniy,t) + jmid{y,t) + jfi{y,t), dSOj), dH]), dSSj) and dM])- 
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Proposition 2.19. For t > 1, 


(74a) 


(74b) 

// N 1 kA 

jm^d{y,t) =0{1) I (a ^(1-a)* ^ + '^{1- afk^ ^ + (1 - 77(t) 

\ i=i fe=i / 

k=l ^<s<t 

(74c) 

j/z( 2 /,i) =0(1) sup |j(y,t) - j( 2 /',i')l+0(1)(1 - o)^"’ {j{t)+ \\9in\\^A 

y.y^dD 


Proof. 

For jin, since ||| < ^ 


|jm(y,OI 


'l€l< 




/ O'TT \ ^ 


dxd^ 


= 0 ( 1 ) / 
-'iSKf 


77(t)M(0 + ||5*n|lc 


Vi (l + KI)'^ 


dry 




0 ( 1 ) 


^77(t) + IlSmlloo,;^) • 


By (112]), 


1 A r fe 

jrmd{y,t) < |a^(l - --ti)] 


(i^) ^ ^-(x«)(0 + (i-)'= (^aE(i-a)-^7(y,0 ^ 




VO 
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Direct computations yield 


"' tiL'-' ^ ' *w.,,(o 


-9in{x(k) - ... -tk),^’') 


K 


= 0 (l)^(l-a)^ 


fc=l 




hi' 


loo,// 


dxd$, 


( K 


0{l)(j{t) + Un\h,, 




K 


E(l-«)"^ ifd = 2 . 


tfc=l 


It ia easy to show that 


t-h- ...-U>-t 


k > 2i on A 


k-) 


and therefore, we rewrite I as 
K k 


fc=l' 


2=1 


27r 




^i?r(x(,))^ 

K K 

2=1 fe=2 

' LA^/ 2 J 


^22-1 K 


K 


K 


2=1 \fc=2 fc = 22/ 2=[K/2J+1 

Direct computations yield 


k—i 


K 


K 


K 


a E (1 -a)* ^E(--- 

i=\_K/2\-\-l k—i 

K 

= o{\)j{t) E (i-«r'«E 

i=LK/2j-|-l 


k—i 


i ^ 


271 


M{^)d^dy. 


LK/2J 




fc=2i 


E(i-a)* ^«E / (i(y>0 


j<s<t 

= 0 ( 1 ) sup 

i<s<t 


lK/2\ 


K 


= 0 ( 1 ) sup (llillioo) (s) E (!-«)* ^“E 


i=l 

1 ^’ 


k—2i 


— 2o ^ ^k 


271 


27r 


i?T(x(,)) 


M{^)d^dy 


^T(x(i))(l)d€d5* 
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LA/2J 


2i-l 




k—i 


f 27T 


Vi?r(x(i)) 




2i-l 


lK/2j 

= omit) E (i-«r'«E 


27r 


^ \ RT^ 

k—i 




lK/2\ - if d = 1, 

= 0(l):7(t) E 

ifii = 2. 


For jfi, 




/ 271 




log(t+i) 


-{I-a) 5 (x(a:)-| (t - ti - ... - tx), $, ) jdxd| 

= o(i)(i-a)^(:rw + ||g»IL,J, 


27r 


l€l>- 


log(t+l) 


. -R 

= / KE+E )(■■■) 


1 ^ 1 ^ '^(l)-^( 2 )l 

log(t+l) 


= 0(1) sup |j(y,i) - j(y'.^')l+0 (l)(l-a)^'’, 

i'G(t—-ffr log(t+l),i) ^ ^ 

yyedD 

where 0 < p < 1 is any fixed number. 


□ 


With the estimates (|74a|) . (|74bl) . and (I74cl) . the main task is to study the RHS of 
(ITia . the fluctuation of j. We have treated this in Sections 12.2112.41 and 12.51 The 
fluctuation estimate, the following theorem, follows directly from Corollary 12.131 
on the temporal fluctuation estimate, and Theorem 12 . 141 on the spacial fluctuation 
estimate. 

Theorem 2.20 (Fluctuation estimate). Lett' < t. y, j/ G dD, then for sufficiently 
large t'/mn, t'/mN and N, 


(75a) j{y,t) - j{y',t') 

= 0(1) ("^sup (s)"] X y y ^ 




log N\^ ^ f log n 


V n 


(1 - (') 


+ 0(1) (' (-V + „»iV^)log(t- + l) ^ ^ _ 

when d = 2, 
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(75b) 

+ o(i)(^ 


= 0 ( 1 ) 


mrrr log(i' + 1 ) 


i '2 


(t - t') + 


(l-a)"*) {J{t) + \\9^n\\^,,) 

(^) ) ^sup (|j+(s)| + |j-(s)|) 

when d = 1. 


4-<s<t 


Here we will complete the proof Theorem 12.31 using Theorem 12.201 
Recall K rs iog(t+i) as t 1. And hence for each 0 < p < 1, Rp log(t +1) = o{t) 
as t ^ 1. Plugging Theorem 12.201 into (I74c|) . we obtain the following estimate of 
jfi- 


(76a) 

= 0(1) i - ^ -^ + (1-0;) + (1 - a)Vios(‘+i)l j 


+ 0 ( 1 ) 






(log(t + 1 ) 


P 

log(t+l), 


when d = 2, 


(76b) 

+ 0 ( 1 ) 


= 0 ( 1 ) 


m^nr \og{t + 1 ) 


+ (1 - a)^ + (1 - a)(isiif+TT) \ 


((lSi(7TI))'’log(t + l) + t«)\ .^.2 , s 

^ ^^ +(w ( b+(s)l + b-(s)l 

n2 j \t<i/ t/2<s<t^ ^ 


when d = 1 . 


We first establish the nnifrom bonndedness of j: 

Lemma 2.21. The boundary flux j is uniformly bounded: 

(77) 77(t) = 0(l)||g.„|U_^, t> 0 . 

Proof. 

From GM: and dZH), we have for t > 1 , 

0 ( 1 ) 


jm(y,t) =-^ {\\9in\\^,f,+ J{t)) 


jmid{y —^(1) 


(log(l + i))^ 


(ll5™|loo,/i + 77(t)^ 


My. 0 + (1 - «r + (1 - 

I (ll5m|loo,/i + 77(1)^ 


rrdN-‘ log(t+l) ^ _ 2 

0 , d= 1 


+ 


log(t+l) 

*3 


+ ( (T^) + (!=F) ’ ™«) (Eifcj) log(i + 1), d = 2, 

( log(t+l) )^ Iog(t+l)+t'' 




(+■ 


, d = 1 


J{t). 
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Therefore, 

=j^n{y,t) + jrmd{y,t) + jfl{y,t) 

1 


= 0 ( 1 ) 


(log(< + 1 ))^' 


,(■_t_W log(f + 1 ) 

(1 - Q,)lisiIi+T)J +(i_a)™H-- - 


+ ■ 


+ 0 ( 1 ) 


'W 

0 


-,d = 2 , 
, d = 1 


hr. 


+ ■ 


1 +_t_+ log(f + 1 ) 

(log(l+,)). + <l-°>‘'^’ +(l-ar +-i 

'(to) ’ „JV + (!2t!l) (jjto)”l„g(, + 1 ), d = 


log(t+l) 

^(7+T 


So far we only have to assume t/mn, t/Nm, N ^ 1 . Now we set n = n{t) = 
m = m{t) = and N = N{t) = , where 0 < ri,r 2 ,r^ < 1 are to be 

determined. In order to get 


lim 


1 


(1 _ + (1 _ + +++w«(i+i)'|, 


lim 

t—¥00 


we need 


( log(t + l) )^ log(t + l)+t'' 


td +1 

V P 


t-J-oo (log(l + t)Y 

'"’Kllf*" + ( (+=)+«+ (!to)M + ,), j 


= 2 , 


1 “ 

n 2 


(f)' 


d= 1 




t 2 + t 3 +p < -n, 

P<lrs, 

ri+r2 < 1 , 

T2 + T3 < 1, 

r2 < q, 


P,q < 2 ^ 1 - 

This can be done by choosing any 0 < ri <6/7 and setting q = ri/3, r 2 = ri/ 6 , = 

ri/12,p = ri/36. Therefore, there exists t, > 0 such that, for all t > t,. 


(78) 


Il5=) it) ,d = 2 
[\U{t)\l\j.{t)\),d = l 


<0{l)\\g.n\U^ + ljit). 


Moreover, from m, now that t* is fixed, j is bounded by a constant multiple of 
llff™lloo M — t < i*- Hence, (1751) actually holds for all t, and this implies 


that hJ^it) = 0 ( 1 ) Ilmira 11 ^ and the lemma is proved. 


□ 


With the boundedness of j, dZZl), we can perform the second a priori estimate 

to obtain the {at + l)”"^ + (1 — a)*™ decay of j, Theorem 12.31 First, since J = 
O(l)||gm||oo, ^ and the following Lemma l2.221 we can rewrite our previous estimates 
(I74al) , (I74bl) , and (175)) as the following in Proposition 12.231 We will need the 
following identities, whose simple proof is omitted. 


Jit)- 
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Lemma 2.22. For 0 < a; < 1, 


2 ^^^ - (l-x)3- 


fe=l 


Proposition 2.23. For t > 1, 

(79a) jin{y^t) = ll^^y 


Ioo,/^ ’ 


(79b) )„„(„,() =0(1) (^ + (l-°)‘’ (i„g(iV,)). ) Il9.»ll,.„ 


+ 0 ( 1 ) 


(iog(i -h t)Y t 


and 

(79c) 

jfi{y,t) = 0(1) 


[nTn^ + nt'N'^) log(t -I- 1) 
13 


(IIjIIl^) W, 


+ { 1 - a)™ + { 1 - Q()(is5TtFTy) ^ 


/ log TV \ ^ / log 


+ 0 ( 1 ,( s„p^J||;1L,)(4 (^) + 


mN 


f ^ 


\\og{l + t) 

when d = 2, 


log(l+l), 


(179^) 

+ 0 ( 1 ) 


0 ( 1 ) + (1 -„)". + (1 - „)(.=ifai)') b.„iu„ 

((toi+n) log(l + 1 ) + 1 ") /miA / \ 

^ ^-^ + (77 (b+(s)l + |j-(s)l 

n ,2 VI''/ / t/2<s<i V / 


when d = 1. 


Now we are ready to prove Theorem 12.31 
Definition 2.24. The a priori norm of j is a function of t defined as 


(as) ^ + (1 — a)'* 


(80) Af(t) = sup 

0<s<t 


-1 


Jh^is) 


for d=2. 


(as) ^ + (1-a)'*““) (|j+(s)| + |j_(s)|) ford=l, 


where j±(s) = j(±l,s). 

Proof, of Theorem \2.3\ 
From (17^ and for 1 > 1, 


imid{y,t) =0(1) + (1 - a)(i^^ii^) 

+0(1) 


(iog(i + t)Y 


(log(l + 1))'' 

+ {I-aY^\M{t) 


hv. 






















EQUILIBRATING EFFECT OF MAXWELL-TYPE BOUNDARY CONDITION IN HIGHLY RAREFIED GAS 


and 


iMy.O = + f 'Xl 


+ 0 ( 1 ) + (1 - a)‘ 

Therefore, 


M{t) 


( log<l + t) log(l+t)+*‘^ 


+ (f)' 


llfl'inlloo, 

■t)„ d = 2 
, d = 1 


=0(1) -h fl — ^- I 

^ ^ (log(l + t))" 

f t \P ( rn^N^ log(*+l) J_9 
+ (1 — Of)'*°*<*+i)+< ’ 

^ ’ \ 0 , d= 1 

+ 0(1) ({at)-^ + (1 - a)‘™ ) J\f{t) 


1 


^d+ 2 ^d+l 


fd+1 


+ (1 - ay 


\\9t: 


(iog(i + t)y 


(+^) ’ + (i+) ’ X (Eifer)" log(l + i), d = 2 


( log(*l+t) )^ l°g(l+*)+*‘* \ , /'m\2 

-+ 


, d = 1 


Now we choose 0 < ri,r 2 ,r 3 <C 1 and set n = n{t) = m = m{t) = 

N = N(t) = so that 


(1 - a)( 


log(t+l) , 


m‘'+ 2 ^d+i log(t + 1 ) / _ t _ 

- ld+1 - - + (1 “ “) + (1 - 

= o((at)-‘^ + (l-a)*^), 


(log(l + t))^ 

log(t+l) ^ _ 2 

^ ’d -1 
^2 , a — 1 

’(¥)* + (+=)*”>«) (Eifer)" >“g('+ *)>‘‘ = 2 

(iog(i+t) y iog(i+*)+*' 


lim --^ 

t-yoo (log(l + t)y 


riy 


= 0 . 


(sr 


d=l 


This can be done if 0 < ri,r 2 ,r 3 ,p, and q satisfy 

1 


T2 +T3 +p < -ri, 

P<ir3, 


< r 2 , 


ri +r 2 < 


T 2 + T3 < 

r 2 < q, 

1 


1 


d + 2’ 
1 

d + 1’ 


p,q < 2 ^ 1 - 

In fact, we may choose p = T 3 = ^,T 2 = ■^,q = ^,ri = jh so that the above 
inequalities hold. Consequently, we can find some sufficiently large t* > 0 such that 
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for all t>tf 
(81) 

(IIjIIu-) W ,d=2 


(|j+(i)l + d 


= 1 


< ((af)-'^+(l-a)‘™ ) ( 0 ( 1 ) ) . 


= 2 
= 1 


Therefore, 

(82) 

((af)-2 + (l-a)‘™^”'(||j||^^)(f) ,d 

+ (1 _ af^y ( \Mt)\ + |j_(f)| ), d 

Moreover, by Lemma [2.211 for all t < 

(83) ^ 

((af)-2 + (l-a)‘^^”'(|lj||^^)(f) ,d 


< 0(1) llff™|loo,^+2-^(^)- 


= 2 


t) i + (l-a)‘™) (|j+(f)| + |j_(f)|), d= 1 


^ =0(l)||ff„|| 


00,/2 ' 


Hence, (1^ actually holds for all t, and this implies Af{t) = 0(1) Hffmlloo fj,- From 
this estimate, the definition of A/’(t), (1501) . and the boundedness of j{y,t), (1771) . it 
is easy to see that 


j{y,t) = 0(1) ||5„|1^ 


1 


\ {l + aty 


+ ii-ay 


and the theorem is proved. 


□ 


3. Damped Free Molecular Flow 

In this section we consider the damped free molecular flow, with a goal of proving 
our main theorem. Theorem 11.51 and Theorem ll.61 As noted in Section 1, we treat 
nondimensionalized Boltzmann equation so from now on we set T* = 1. 

We first review some basic properties of the collision operator Q and the lin¬ 
earized collision operator L. Both operators act on ^ but not on x. Hence we will 
frequently neglect the spacial dependence in the following discussion. 


3.1. Preliminaries. 

Recall the collision operator Q(-, •): 


Qig, h)iC) = \ j {giOHO + KOgiO - g{C)KU - KOgiu) 

S2xR3 

xB( 0 ,|C-CI)rfMC. 


As intermolecular collision conserves total molecular number, total momentum, and 
total energy, we have the following identity, cf. m. 


(84) 



Q{g,hm dC = 



Ij Cl ICP are called the collision invariances. 

In this paper we assume an inverse power hard potential with Grad’s angular 

11 — 4 

cut-off or hard spheres. Under this model, B{9, |C* — Cl) 1C “ C*l~^l cos0|, for 
some M > 4. For the linearized collision operator Lf = As a 
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Vm(C)(L/)(C) dC = 


direct consequence of 

(85) / ( ^ ^ 

|^|2 

Recall that L can be decomposed as the difference of an integral operator K and 
a multiplicative operator v: L = K — v. Moreover, we have the following estimates, 
cf. 0 : 



( 86 ) ^( 0 ^( 1 + |C|)^ 

(87) 7f(C,CJ=0(l)IC-C 

And it is well-known that 
( 88 a) 

( 88 b) 


^-4 

u e 8 


0 ( 1 ) 




K-C 




Q 


v^fM 




y/M 


=0(1) X IIV'IIl- + X 


Let 


= esssup(l + Id) 1/(01- 

« C 6 R^ 


Also, for any 7 > 0, 


(89) 


\\Kf\\ 


T 00, —-y — 


IfiUKiCXMC 


= 0 ( 1 ) 2 '> 




(91) 


K 


(// 


r 00,— r 


, for e <C 1 . 


(1 + id)^ 

(IRiH) with 7 = 0 states that K : —>■ is a bounded operator. Actually, more 

is true: K : —>■ is a bounded operator for each ,5 > 0, 0. However, 

(IRiHl suffices for our purpose. From (IMll . 

(90) II^/IIl-= o(i)||u/||^^.. 

Let feiC) = l{|c|<e}(l + ICI)’^ for e < 1 . By ([57]) we have 

0(e) , d = 1 

0(e2|loge|), d = 2 

The details can be found in m- 

As mentioned in Section 1, we take (HU as our linearized problem. To show 
an exponential decay property of (HU, we conduct two reductions: first we reduce 
(dt + J2Cidxi - ^L) to (dt + + ^), and then reduce {dt + + ^) 

to {dt + YO^Xi)- To facilitate the following discussion, we invoke the notion of 
solution operator. 

Definition 3.1. Sf^ are linear operators defined as the following: 

(92) (sLB(/„))(x,C) = /(x,C,t), 

- -Lf = 0, /(x, C, 0) = /„(x, C) 

OXj K, 


Maxwell-type boundary condition HOD 
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(93) (s°^(/„))(x,C) =/(x,C,i). 

|f + g4^ + r-^ = ». /(x,<,0) = /.„(x.C) 

I Maxwell-type boundary condition m 


(sp (gm))(x, C) = ff(x,C,i)> 


^ + = 0, gi^,C, 0) = 5m(x,C) 

i—1 

Maxwell-type boundary condition 


We call and Sp^ the solution operators for the Linearized Boltzmann 

equation, free molecular flow with damping, and free molecular flow, respectively. 


From ®, any /r > 4 is an admissible choice for free molecular flow. Since there 
is no need to vary for definiteness, from now on we fix ^ = 5. The pointwise 
results for the free molecular flow, Theorem ll.il is written in the following form for 

^ “ 400 ■ 

Theorem 3.2 (Main Theorem of Free Molecular Flow: Solution Operator Form). 
For fin e 0 < 7 < 1, 


Sfr (uf.nVM) (xX) 

Vm{i + \cx 

Sp- (/»yM) {xX) 

VMii + \c\)'y 


=0(l)ll/™|loo,-^KC), 

=o(i)||/™IL,_. 


If, in addition, f fin\/MdxdC = 0, 


Sp(/»yM) [xX) 
v^(l + |CI)^ 

= 0(1) WfinWao,--/ I (jY+Xdj^ + 

So far we have only obtained the existence, uniqueness, and pointwise esitmate 
of free molecular flow. Theorem o Now we will settle down this issue for the 
damped free molecular flow, (1931) , and obtain a pointwise estimate of the solution. 


- a) 




If , 1 +]ir 

L 1400 J L 


l^l< i 

t 400 J 


3.2. Global Existence and Boundedness. 

In this subsection we establish the global existence and boundedness of boundary 
flux of the damped free molecular flow: 


(94) 


+ - 0 . 

2 — 1 

g‘'(x,C,o) =5P„(x,c) 

Maxwell-type boundary condition (|51). 


This can be done easily by the comparison with the free molecular flow. 
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By the characteristic method, solutions of (IMl) can be represented as 
(95) 


/ m — 1 




^(l-a)"e (x(.+i),t - h - it2) ( ^r(xo+i)) 

+ (1 - 1)^2), C") for h < t, 


<'«) 


. e - 5*n(x-|t,C) fort<ti, 


where 


Consequently, 


lilt - |X-X(1)| 

m = [—-+ 1 . 

|X(1) -X(2)| 


fiy,t)= j (-li-n(y))e Ci)rfCi 


t< 


|y-y{i) I 
l«il 


(-I1 ■ n(y))ae“ 


.idiii. 


t> 


|y-y(i)l 

kil 




t> 


|y-y(i) 

ISil 


(-€1 •n(y))(l-a)e ] dC^. 


One can follow the discussion in Subsection 12.21 to derive the formula for 


n / k k ^ 

^) = H f{k) (y. + Y. HkY-M) (y’ 

/c—0 \ 1—1 ki-\-...-\-ki—l y 

n+1 n+1 


+ '^ra+l(y)0 + E E 


i/'(Z,n+l) 

{ki,...,ki) 




1—1 ki-\-...-\-ki—l 


where 


j(k)iy^t) = 


_ •'(ci) ly(i)-y(i-n)l 

e *=“ 


l «1 


„ ‘Ui |y(o-y(i+i)l iy(fc)-x(fc+i)i 

^ I>;i I ^ \ck I 


l^il 


(-I1 ■ n(y)) (l-a)'=e 


i«fi 

|y(0-y(i-n)l 

=0 


14^1 


^27) 




2=0 


irii 
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= / n • n(y(fci,....fc._i)) (1 - a)'"' 

. n i—1 


A(‘.k) 

(ki . k,) 

_ ^ ly(fcl.fci-1.3-l)-y(fcl.fci-l,3)l 

e ^=1 " 


27T 


{1-af 


-ki-...-ki 


k-ki^--ki |y(fcj,,,,,fc,,i_i)-y(fcj,,,,,fc,,i)l^ 

-S--) 


ISLi I 


... ^ ^ .fe._^_,_i)-y(fc^.fei_i,j)l '‘My(fci.....fc,,i_i)-y(fci,...,fc,,i)l 

- ^ ^-^ - 

z=lj=l I 1 = 1 l^z + l I 


1/ I ^k—k\ — ... — ki(_L ly(fci,...ifc i-iJ —1) y{ki,...^ki—i,j)\ 

9in I y(fci ,...,fei,fc—fci—... —fci) ^l+l / , / , 

i=l j=l 


le 


i-i| 


"■ |y(fci.... ,ki,i-l) y (fci,...,fci,z) I ^k-ki-...-ki\ 7* 7> 

1| /’>/+! I ^Z+1 * • * >1? 

i=i l^;+il / 


_ " .7(g«) |ym-y(i+i)l 

U%i(y,0= J (-^i-n(y))e -o » (i _«)"+! 


‘>E 


|y(i)-y(i+i)i 


17 / f _ V y(»+i)l >n+l\ ,> 

? ly(n+l))^ / 7 l^i I >>1 I 


.17(Z.«+1) 7 7 _ 


L 

Y[{-^^■^iy(k„...M-l)) (1 


— a)*‘ 


:,(Z,n4-l) 


h •'«{ . ki_i,j-l)-y(ki .fei_i,3)l 

“ ^ —;;- 771^. - 


my(kA....M))) ^ 


^n+l '“1 ■■ '^l |y(fc^.fc,,i-l)-y(fci,...,fc,,t)l 

\\ " |y(fci,...,fc i-i j-i) ~ y(fci,...,fci_i,j) I 

I y(fci,...,fci)>^ / 7 / 7 


i=i j=i 


le 


j-ii 


n-|-l —/ci —... —fc; I I '' 

E iy(fci,...,fci,i-i) ~ y(fci,...,fci,j)I >n-^l —fci —... —fc; 

I /’i —11 ’ >Z -|-1 

i=l 


ICi 


^Ci+i • • • dCi, 


and 


y(o) — y? y{ki,...,ki,0) — y{ki,...,ki)'> 

y(ki,...,ki-l,i) ~yB (y{ki,...,ki-i,i-l)i l^pd) ’ 

Cl = Cr^ - 2(ir^ • n(y(fc^_,„_fc,_^_,)))n(y(fc^_ 
Ci = {C,rii), 
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A 


(i,k) 


I ki 






i-l I 


a k - ki - ... - ki = 0, 


< 


i=i j=i 

| y ( fci ,..., fci ) ~ y ( fci ,..., fci , i ) I 
l^;+i I 

\y (ki^... ,ki ,k—ki —.. . — ki) y {ki. .,ki ^ki ,k—ki —.. . — ki + l)\ 


2=1 




ii k - ki - ... - ki > 0 


B 


(Z,ra+1) 

(ki,...,ki) 

I ki 


^ ^ y^ y^ |y(fci,---,fci-i,j-i) y(fci,...,fci_i,j)l ^ ^ |y(fci,.. .,fci,i-i) y{ki,...,ki,i)\ 


i=l 3 = 1 




i=l 


le 


z-l-l I 


Since the kernel of free molecular flow always dominates that of damped free 
molecular flow: 


j'(k)iy^i) < 


(-€i • n(y)) (1 - a)* 


|y{i)-y(i+i)l ^ |y{fe)-y{fe-,-i) I 
^<^- 2 = -fTO-<-TTJE- 


isil i«fi 

* - d(( - E d'l 


2=0 


111 I 






"ny,...„)(d'"') (-«!+. ■ «(y(...h))) (1 


— fei ——/cj 


I 'y{k-i,...,ki,k—k\— ... — ki) 


t k' I 

k-ki-...-ki _ y^ y^ |y(/ci,...,fc i_i,j-i) ~ yj') 

i=l 3 = 1 


14^ 


i-l I 


k-ki-...-ki I _ 

|y(fci,...,/i;i,2—1) y(fci,..., 

2 = 1 


2-l| 


l^^ + l 


rfCz+i-.-dCi- 


oo / k k 

i"(y’ = I] iw (y’ 0 + 51 51 (y^ 0 

k—O \ /=! ki-\-...-\-ki—l 

= o { i ) MJU ,. 

This proves global existence of the solution and uniform boundedness of j'^ by 
comparison method, for the case > 0. For general initial configuration, let 
g'^’^ be the solution of damped free molecular flow, (IMl) . with initial configuration 
the positive/negative part of Let be the flux of g"’^. Note that 
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(IMl) is linear, j'^ = ■ (However, ^ in general.) 

ir(y,i)l 

=o(i)(||{<?rj+L,5 + IKff»)”L.5) =o(i)iiffrniioo.5- 

To sum up, we have 

Theorem 3.3. For all t ^^6 solution of (IMll exists globally, with 

{\\r\\Lr)it)=Oil)Kn\L,,- 

Therefore, for fin G 

l|SP''''(/.n)L,_^ = 0(l)||/zn|L,_^. 

3.3. A Pointwise Estimate. 

For convenient, from now on we will frequently abbreviate functions /(x, ^,t), 
ip{x,<^,t), etc., as f{t), t{j{t), etc.. 

Lemma 3.4. Suppose that fin G , for some eonstant 7, 0 < 7 < 1. Then, 

under the zero total initial molecular number assumption J fin'/Mdxd^ = 0, 


(1 + KI)'^ 




( — I — 

V(1 + a<)^ 


+ (1 


Proof. 

By Duhamel principle, (1931) is equivalent to 




3yy 


} 



399 I 
tiM J 



/(t)=SP^(/„) = 

From Theorem 13.21 


Sf^if^nVM) 1 f* Sfl,iufis)VM) 


'J M K Jq yf M 


ds. 


aFr 


(/« 


v^(l + |C|)^ 


= o(i)ll/»IL.-. 


1 . . tTOCT , ^ 

+ (1 - a) 2 II 


(1 + at)'^ 


I l$l>-^| 

L t400 J 


Il€l<-w| I 

I. t 400 J '' 


From Theorem 13.21 and Theorem 13.31 

SfL,ii^f{s)VM) 


Vm{i + \<:\p 


= 0 (l)||/(s)|U.-^u(C) = O(l)||/,„||oo.-^u(c). 


Hence 


1 /■* sp:,(u/(s)^) „ 

Wo yM(i + |cl)'^ -o(i)^||/™||oo,-WC). 


□ 


Our next step is to remove the undesirable factor of u(^) in Lemma Td.41 To do 
this, we conduct a posteriori estimate through the characteristic method. 
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Theorem 3.5. Suppose that fin G 0 < 7 < 1. Then under the zero total 

initial molecular number assumption f fin\/Mdxd^ = 0, 


SP^(/.n) 

(1 + KIP 


0 (l)||/.n|L.-. 


/ 1 / 1400 / 1 


xlr 1 

+ 1 f 1 

t \ 

+ - 

\ I€I>-355- ( 

L 4400 J 

■ \ ICI<-3gg- ( 

l. t400 J 

> K 1 


Proof. 

By the characteristic method, 


sP^''(/„) (x,C) = /(x,Ct) 

' m —1 

a ^(1 - a)*e (^X(^._^i),t - ti - 

i^O 

+ (1 - a)"^e ^‘/m(x(m) - -h- {m 


it2) 


27T 


i?T(x(,+i)) 




-l)t 2 ),C™) for ti < t, 


From Lemma 1311 

f 1 s 4M / 1 

j(y..) = 0 (i)|l/,„IU_,(^i^-^ + (i-a) . +(^ 

Hence Theorem id. 51 follows. □ 



4. Steady State Solution of The Boltzmann Equation and Its Time 

Asymptotic Stability 


4.1. Linearized Boltzmann Equation. 

In this subsection, we study the linearized Boltzmann equation (HU. We first 
establish the local (in time) existence and a local estimate. By local in time we 
mean L not t <C 1. From now on we always consider perturbations of the 

form /Vm {fin G 

Theorem 4.1 (Local Existence and Estimate). Let fm G L'^^~~^ , 0 < 7 < 1, then 
there exists a constant c* > 0 , such that whenever - < c*, the solution of (HU) 
exists and satisfies 


Sr(f^n) 
(1 + KI)^ 


= o(i)|l/„|l^,_. 


f ^ 

[(l + at)'^ 


t 400 f 1 

+ (1 - a) 2 -I- -- —— 

V (1 + t)*oo 

Xlr- -i-l-llf 1 

\ l^l> 399 f { l€l< 3SS f 

1. tlM J L tiro J 



Proof. 

Write (ITTl) as 

(96) f{t) = SL^(/„) = SP^(/„) + - fsfJf{Kf{s))ds. 

Jo 

We solve (IMll by iteration: 

/(o)(t) ^ SP^(/„), /«(t) = SP^(/,„) + [Kf^^-^\s)) ds, i > 1. 
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From Theorem 13.51 we have, 
(97) 


(1 + KI)" 


= o{i)\\u\_ 


00,-7 


(1 + at)‘^ 


1 

+ (l-a) = 

4 


\ d-\-l 

, (1 + t) 400 

J 

X 1 

\ l€l>-^ f 

4 400 J 

+ 11^ 

[ l€l<-w| 

L 4400 J 

4j 


From (1^ . 


= 0 (i)||/„|l (1 + -). 

00,-7 \ 


Consequently, from (IMl) and Theorem 13.31 




00,-7 N Jq 


jDFr 

^t—s 


(iF/(°)(s)) 


ds 


00,-7 


0 ( 1 ) 


I (i+0 (ow 


Similarly, we obtain, by induction, 


=||/,„|| (i + :i) ( 0 ( 1 )-) ,*>!,. 


00,-7 


iioo,—7 


Hence, for ^ sufficiently small. 


^||/W(i) =0(l)-<oo, 


and Theorem 14.II follows. 


□ 


With the local estimate of Theorem HU we are ready to prove the global (in 
time) exponential decay of dTlT) . Theorem 11.31 Recall that vq = inf u(^). 

Proof, of Theorem \l.S\ 

For simplicity we write S^®(/i„)(x, C) as f{t) in this proof. Define 


P{t) = sup e - ||/(s) 

0<s<t 


loo ,—7 ' 


Let c < C:4:/2 be a small constant to be specified later. Recall that c* is a constant 
given in Theorem 14.1[ Since 2c < c^, we may apply Theorem [TT] to obtain 

Pi2cK) = sup I — : X e D, C e 0 < t < 2 ck| 

< sup ||/(t)|U sup ||/(i)lloo ,-7 = ^(1) ||/m|L,-,y 

0<i<2cK 

By the definition of P, it is clear that P increases with t. We now claim that under 
some appropriate choice of c and k 


(98) 


P{t) < P{t — ck), whenever t > 2ck. 


00 ,- 7 ’ 


This claim implies P{t) = 0(1) ||/m|| 


which proves this theorem. 
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Tracing back from time t to the earlier time t — cnhy the characteristic method, 
we can represent /(t) as: 

(99) 

/(x, C, t) - c<, C, t - ck) 

1 .(O 

H— / e K{f{s))(x — s^,C,t — s)ds, when > ck, 

^ Jo 

1 t / 9 ^ 

/(x,C,i) = X! 


2 = 0 


i?T(x(,+i)) 


^T(X(,+ 1 ))(C) 


+ (l-a)*- 


.KO 




e - ‘'K{f{s)){yi(^i) - sCX\t-s)ds 


+(1 — a)"^e — ti — (m — 1 )^ 2 ), C™) t — ck), when ti < ck, 


where 


ti —■ 

t2 =■ 

m =[ 


^( 1)1 


III ’ 

^( 1 ) -^( 2)1 

III ’ 

|||CK- |x-X(i) 

|X(1) -X(2)| 


J +1. 


Consider first j(x(i_|_i),t — ti — 1 ^ 2 ) with ti + it 2 < ck, for i = 0,...,m — 1. 
We trace back an extra 2ck — ti — it 2 amount of time to arrive at t — 2ck. Since 
ti + it 2 < CK for i = 0,..., m — 1, by Theorem 14.II we have 


( 100 ) 

j(x(,-Hi),t-ti-it2) =0(1) / -I* ■n||/(t-2cK)||^ _^(l + |C,|)^ 

J •n<0 


y'^{2ck — ti— it2)’^ 


(2c^-t, -it,) 400 

+ (1- a) 2 + 


(2ck — ti — it2) 


X 1 




l€J> 


-'-!UU 

(2cK-ti-it2)OT 


I^J<- 


( 2 cK-ti - 2 * 2 ) 400 


IM" f J 


= 0 ( 1 ) ||/(t- 2 ck)| 1 ^^ . 


{ack) 


1 GidiSk 

+ (1 — a) 2 -p 


/ \ 399 

(CKj 400 


d +1 


+ c . 



Therefore, 


m —1 


^(l-a)T 


-^C^l+*^ 2 ) , 


2=0 


27r 


1 - i, - it,) 


= 0(l)||/(t-2cK)|L,_, 


{ack) 


1 . (ctt)WU 

+ (1 - a) 2 + 


/ \ ^9 

(CK,) 400 


d +1 


27r 

mi 


M{C). 


For Theorem 14.II to apply, we need f f(t — 2cK)'/MdxdC = 0. This is true because 
of (( 551 ) . 
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Consider next e '^^K{f{t — s))ds. We also trace back to the time 

t — 2ck. From Theorem I4JI (15^ . and (lOTl) . 

( 101 ) 


E (!-«)*- 




(O 


2=0 




e - "K{f{s)){x^i)-sC,C,t-s)ds 


=0{l)\\f{t-2c^)\\^_^{l + \C\r 

""- 1 I |•tl+...+ti+l ( ^ 


E- 


2=0 


Kytj+...+t. la'^( 2 cN-s)' 


(2CK-S) 400 

+ (l-a)^---f 


1 


/O \ 

[Zen — s) 400 


d+l 


+ 


1 


/o \ 

[2cn — s) 400 
2 


, when d = 1 


/n \ 

[Zen — S) 400 


log{2cn — s), when d = 2 


1 


> e^ds 


=0(1) ||/(t - 2c«)|l^,_^ (1 + Id)" -P(cn) + , 


where 


P{z) = 


1 zWU 1 

- \-z[l — a) 2 +2:400 for d = 1, 

a 

1 . . z At log z 

^ — h z(l - a) 2 -I-for d = 2. 


199 
2 : 200 


Plugging (llOOp and (jlOll) back to (|5^ , we have 

(1 + KI)" 

g-(^o-<2l)cjr(t - ck) + C'e^^<T{t - 2 ck) f-P( ck) + 


for ti > CK, 


< 


C'e^‘=‘'^P(t-2cK)< 


(ack) 


1 ^ 

+ (1 — a) 2 + 


/ \ 
en) 400 


d-hl 


+ I + (1 - - ck) 


<P{t — ck) 


X 




1 - a + ■ 


for fi < CK, 


for ti > CK, 


(ack) 


1 . (cK)m 

+ (1 — a) 2 -f 


. \ 399 

[en] 400 


d+i 


+ c 


+ ( —P{ck) + c‘ 


for ti < CK, 


for some positive constant C'. 

For (IMl) to hold, we need 

( 102 ) 




l-a + 


1 *°° 

+ (1 — a) 2 -p 


(ack) 


/ \ 3^ 

[en) 400 


d+l 


+ c + -P{ck) + & 

K 


> < 1 . 
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For this purpose, large k and small c are desirable. We will fix c according to v[ 
and a, and find the admissible choice of k after c has been specified. We fix some 
small c, so small that 


(103) 


c < a, 

i min{uo — u(,a} . 


It is not difficult to see that we can choose some c ~ a(z/o — i^'i) to meet all these 
requirements. From (11031) . 

\K J 4: Hi 




(ack) 


1 (CK) 400 

+ (1 - a) 2 + 


/ N 399 
CAu 400 


d +1 


+ c + ( -P{ck) + P 


< 1 - - + 


1 


(ck) 400 

+ (l-a)^^ + 


{ackY 

For this specific c ^ Q!(r'o ~ need k to satisfy 


/ \ 3^ 

(CK,) 400 


d +1 


H— P{ck) 

K 


(104) 


—P(ck) < c——— ^ a(r'o — u()^, 
K 4 




1 ■ (e,.)TOT 

+ (1 — a) 2 -I- 


(ack) 


1 


/ \ 399 

(CK) 400 


d +1 


H— P{ck) 

K 


a 

< — . 

- 2 


It is not difficult to see that there exists a positive constant Ci, independent of 
c, such that (IT^ implies (11041) . Consequently, under the assumption (fT^ . (IMl) 
holds. □ 


4.2. Exponential Convergence for Pull Boltzmann Equation. 

Using the exponential decay of Linearized Boltzmann equation. Theorem 11.31 
we are able to establish the existence of steady state solution and the exponential 
decay for full Boltzmann equation, © or equivalently dH]). In terms of (0 is 
equivalent to 

(105) fit) = SL^(/„) + i (l (^)) ds 

Y gLB ( Q{S-M + VMf,S-M + VMf) 

K Jo \ VM 

In view of (I105L the equation for the steady state solution $ = <l>(x, ^) is 

^ [ Vm 

We first prove the existence of the steady state solution $, Theorem 11.51 
Proof, of the Theorem 11.51 
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We solve (11061) by Picard iteration: 




jLB 


Q{S-M,S-M) 


ds, 


(107) 

'gLB ( Q{s-M + S' -M + ' 


ir 

K do 


\/M 


ds. 


From ((84ll and (1^ . f ^OigJxdC = 0 for all i. Hence (fTKll follows, provided ||$0) ||^ 
converges. 

Claim: Under the assumptions 1 — T, <C 1 and k :$> 1, ||(<i)0) _ $0-i))||^ = 
[0(1)(1 — T,)]*+^, for i > 0. Here we set = 0. 

Note that the claim directly concludes this theorem. 

We prove the claim by induction. First, we follow the proof of Theorem 13 in 
m to conclude 


(108) 




o(i)(i-T,)[o(i)(i-r,)]”+\ 


The next step is to remove the undesirable factor v sitting under (t) on the 
LHS of (110811 . To do this, we conduct a posterior estimate on F^^'>. We note that 

F(») 

can be defined equivalently by the the differential equation 


L/1/ ^ I ^ 

\ *=1 / 


1 


Q + 2(5' - M), 


K y/M 

Maxwell-type boundary condition (nni), 

u(")(x,C,0) = 0. 
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Hence can be represented by the characteristic method as: 


(109) 


X - f (x(,) - rs, C, t-s)ds 

^ Jo 


lf\- 

K Jo 


Q + 2(5 - M), 


^/M 


(x(j) - Cs,C\t- s) ds 


m —1 






2=0 \ 

'ti + ... + U 


27r 








Q + 2(5 - M), 


/M 


(x(j) -Cs,C\t- s) ds . 


From (ll()8L 


m —1 


( 110 ) 


/ Stt \ ^ 

^(l-a)*e-^(‘i+**=)aj(x(,+i),t-ti - Its) 


2=0 
m —1 


=a ^ (1 - a)*[0(l)(l - r,)]"+^ = [0(1)(1 - r,)]"+^; 


2=0 


and, from (jlQSp and 


(111) 


E (!-«)*- 




2=0 




(x(,)-rs,C,i-s) ds 


=[0(i)(i-r,)r"E7 / e-^^ds 

i=0 

= [0(l)(l-r*)r+2i [ e-^^^ds 
^ Jo 

=[ 0 ( i )( i - r *)]"+2 
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From f|88ap . (fT5|) . and the induction hypothesis, 

I e~~^ 

Hi 




Q + 2(5' - M), 


^/M 


ds 


( 112 ) < 


— e 


Q + 2(5 - M), 




=0(1) +o(i)(i-r,)) 

= [0(1)(1-T,)][0(l)(l-T,)r+1. 

Plugging (IllOp . (imp , and (11121) back to (11091) . we obtain 

U”)(t) = [0(l)(l-T,)]"+2. 

Consequently, = [0(1)(1 - r,)]"+2. This concludes the claim 


and therefore this theorem. 


□ 


Now we have already obtained the steady state solution Foe, = 5 + \fM^ for 
full Boltzmann equation ©■ To establish the nonlinear stability for the initial¬ 
boundary value problem (HU), Theorem 11.61 one can follow the proof of Theorem 
14 in m- The idea is basically the same and we omit the details here. 
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